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SUMMARY 


The present report gives the results of a detailed study of the 
flutter characteristics of four representative aircraft wings . This 
study was made using the electric analog computer at the California 
Institute of Technology. During the covers e of this investigation eight 
important parameters of each wing were varied and, in addition, the 
effects of mass, inertia, pitching spring, and location of a concentrated 
mass were investigated for all four wings and several sweepback angles. 

Hie Introduction of this report discusses in general terms the 
flutter characteristics of airplanes. The second section contains a 
discussion of the electric-analog principles that made a study of thin 
magnitude feasible. Hie third section contains a discussion of the aero- 
dynamic and structural approximations made for simplifying the flutter 
analysis of a wing. Hie fourth section gives information relating to 
the errors Introduced by the finite-difference approximations to continu- 
ous aeroelastlc systems. In addition, data are given pertaining to the 
flutter characteristics of a swept-wlng wind-tunnel model and the results 
of computations based on two assumptions regarding aerodynamic forces on 
a swept wing. Hie fifth section lists all pertinent data relating to the 
four representative aircraft wings and the sixth section contains the 
computed flutter characteristics of the four wings. 


introduction 


Flutter is a phenomenon which is observed in the transient or 
unforced response of an aerodynamic system. Mathematically speaking, 
it is observed in the solution of the homogeneous differential equation 
describing the behavior of an airplane in flight through still nonturbulent 
air. An airplane wing which is considered to be a continuous beamlike 
or platelike structure has an infinite number of degrees of freedom, and 
the characteristic equation which describes the transient response has 
an Infinite number of roots. Experience has shown that only the roots of 
lower magnitude (frequency) exhibit the problem of instability or flutter. 
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It Is this fact which mates It possible to predict flatter using an analog 
computer "which represents only the lower frequency modes of the structure 
or using a few normal modes in either digital or analog computation. 

The exponents In the transient response of a linear system are the 
roots of the characteristic equation. Since the characteristic equation 
Involves real parameters , the roots are real or occur as complex conjugate 
pairs. The latter roots are the ones of Interest here. The real part of 
a conjugate pair is the reciprocal of the time constant in the transient 
response and the (positive) Imaginary part is the frequency of oscillation. 
This is illustrated in figure 1. Mathematical description of the tran- 
sient term is 

y = A ie (« to > t + A 2 e (<r ' 1 “ )t 
or in terms of real functions 

y 1 = Ae°^cos(ajt + 0) 

If the real part of the pair of roots a is negative the "transient" 
dies out and the root is said to he stable. If the real part is positive 
the transient grows exponentially until limited by nonllnearitles or 
destruction, and the root is said to flutter. The terminology is not 
strictly correct, hut it is common practice to refer to the exponents of 
the transient response as flutter roots, since they are numerically equal 
to the roots of the characteristic equation. Throughout this report such 
terminology will he used. 

Damping of flutter roots may he measured by two dimensionless num- 
bers £ and g, which differ from each other by a factor of 2. The 
farmer is generally used by control-system engineers; the latter, by 
flutter analysts. They can be defined by the equation for the particular 
term in the transient response given earlier 


y *» Ae°^cos(ajt + 0) =* Ae cos 

SDn-fc 


«= Ae c cos4 


]j(l - C 2 )<%t + frfj 

- («7 4 )] 


ttfat + 


Flutter computations are usually centered around regions where the v alue 

of g lies in the range -0.2 < g < 0.2. In such cases the factor \Jl - 
differs from unity by less than 0.3 percent. For this reason it is custom- 
ary to omit this factor in the trigonometric term giving the following 
approximation: 
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- g ■ 

y « Ae • cos^ooot + ^ ) 

This practice vill he followed In this report. For damping which is 
small, an approximate irule of thumb Is that the damping factor g Is 
nearly equal to the per unit decrement per cycle divided by if. If per- 
cent decrement per cycle 5 Is used, there results the convenient 
approximation 


lOOit 

The flutter roots of an airplane are complex functions of all 
geometrical, structural, and Inertial properties of the airframe as well 
as of the airspeed and air density. With all other properties held 
constant, the lowest airspeed at which the flutter root exhibits neutral 
stability Is called the flutter speed. If g is plotted as a function 
of velocity, the abscissa (speed) at which the curve first crosses the 
axis g = 0 Is the flutter speed. In this study such curves were used 
to determine the flutter speed, but such curves are used In this report 
only to illustrate the behavior of some unusual flutter roots. A tabu- 
lation of flutter speeds does not always give a good picture of the flutter 
characteristics. An example is shown In figure 2 , where the damping of 
two roots Is shown. One root becomes unstable at a speed of about 
300 miles per hour and the other, at a speed of about 600 miles per hour. 

If a parameter variation Increases the damping g of both roots by 0 . 03 , 
one flutter speed Is raised to 350 miles per hour, a 17-percent Increase, 
and the other Is raised to 603 miles per hour, a 0.3-percent Increase. 

A further Increase In g of 0.02 will raise the second flutter speed 
O.lf- percent, to 603 miles per hour, while the first root will now exhibit 
no flutter. It should be emphasized that even though a design speed of, 
say, 300 miles per hour has been surpassed, the system may still be 
regarded as unsatisfactory. A system so close to flutter at a speed of 
360 miles per hour might actually flutter because of weight (fuel) vari- 
ations or minor differences In stiffness resulting from variations within 
the manufacturing tolerances. From the standpoint of this report, an 
three of the sets of roots discussed above will be regarded as having 
essentially the same "flutter characteristics," even though they exhibit 
radically different theoretical flutter speeds. Enphasls Is given to this 
point because remarks to be made later In this report may be misunderstood 
without a clear conception of this viewpoint. 

This Investigation was conducted at the California Institute of 
Technology under the sponsorship and with the financial assistance of the 
National Advisory Committee for Aeronautics. 
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SYMBOIfi 


constants 

half chard 
half chord at root 

half Chord at tip 

symbolic representation of circulatory component of lift force 
due to angle of attack 

lift coefficient 

wing station from root. In. 

Young's modulus of elasticity 

equivalent beam flexural rigidity, (lb) (sq In.) 

experimental j used as a subscript 

flutter frequency, cps 

normal mode frequency of cantilevered engine and nacelle, cps 
flutter frequency for continuous structure 
shear modulus 

equivalent beam torsional rigidity, (lb) (sq In.) 

8%* 

* O 

damping factor of a damped sinusoid, e cos cot 
vertical deflection, positive down. In. 

moment of Inertia per unit length, lb-sec^ < 

torsional stiffness 

Increase In stiffness, percent 

radius of gyration, in. 

semispan of wing 
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twisting moment about elastic aaela per unit length of wing, 
positive nose up, lb 

mass per unit length, lb-sec 2 /sq. In. 

mass of concentrated mass 

fuselage mass 

total wing mass, Ib-sec'Vin. 
total wing mass outside of fuselage 


m lumped mass 

P * P 1> P 2 • Lift force per unit length of wing, positive nose down, lb /in. 
p Laplace transformation variable 

qLn dynamic pressure based on normal component of velocity, 

(l/2)pv n 2 , Ib/sq. In. 

t time, sec 

v alrstream velocity, In. /sec 

▼b flutter velocity of airplane with bare wing 

Vf alrstream velocity at which flutter occurs. In. /sec 

v n com ponen t of alrstream velocity perpendicular to elastic axis, 

v cos A, In. /sec 

v Q reference velocity, in. /sec 

Yqj flutter velocity for continuous wing 

v distance measured along wing 

Xq distance from midchord aft to elastic axis. In. 

distance from quarter chord aft to elastic axis, Xq + (h/2), in. 

distance from three-fourths chord forward to elastic axis, 
b - x-^. In. 


*2 
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distance from elastic axis aft to center of mass. In. 
general variable 

cell size for finite-difference structure 

a absolute pitch angle about elastic axis, positive nose up, 

radians 

percent decrement per cycle 
per unit critical flmnpi ng 

slope of elastic axis or roll about horizontal axis normal to 
elastic axis, positive tip down, radians 

sweepback angle of elastic axis, deg 

air density, lb-sec^/ln.^ 

wing twisting gradient, da/dtr 

real part of pair of roots 

angular frequency, radians /sec 

undamped natural frequency, radians /sec 


EHSCraUE-ANAIflG METHODS OF FLUTTER ANALYSIS 


The use of electrical analogs for the solution of aeroelastlc 
problems has been discussed in detail in reference 1. The purpose of 
the present section is to summarize the principles briefly. For purposes 
of flutter analysis, the B-fcructural system is assumed to be linear, and 
a linear electrical network is constructed whose electrical behavior 
approximates the dynamic behavior of the linearized structure. For this 
purpose, capacitors are ordinarily used to represent concentrated or 
lumped inertia properties. Inductors are used to represent lumped flex- 
ibility properties, and transformers are used to represent the geometrical 
properties of the structure (refs. 1 and 2). In such electrical analogs, 
voltages throughout the network represent velocities In the structure and 
currents represent forces. Electronic equipment is used to produce cur- 
rents which depend on voltages In the electrical system in the same manner 
In which aerodynamic forces depend upon the velocities of the airfoil. 
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The composite electrical structure can toe regarded as an electrical 
model of the aircraft in the same manner that a wind-tunnel model would 
toe regarded as a structural model. The advantage of this approach lies 
in the relative ease with which one can alter the properties of the model, 
thus performing flutter "computations " with great rapidity. It should toe 
emphasized that the normal modes of the structure are not used as tools 
or elements in the analysis . Qhe analysis consists, in fact, in observing 
the behavior of an electrical model of an aircraft in flight. 

Tha t, behavior which is most readily observed is the transient response 
to a sudden disturbance. This method is therefore similar to the testing 
technique which is ordinarily used for wind-tunnel models. An advantage 
of the electrical method is that tuned pulses may toe used, so that sepa- 
ration of two or more nearly unstable or slightly unstable modes of oscil- 
lation is more readily accomplished. Basic recorded data consist of the 
logarithmic decrement of the response and the frequency of oscillation. 
Flutter speed and frequency for any configuration are ordinarily found toy 
computing the damping g and frequency f for specific values of veloc- 
ity and Interpolating to find the frequency and speed at which g is zero. 


AFFR0X3MA3330HS FOR SIMPLIFYING FLUTTER ANALYSIS 


Structural Representation 

For dynamic analysis of airplane wings of large aspect ratio, it is 
customary to treat the wing as a beamlike structure in tooth vertical 
bending and torsion. It is usually assumed for simplicity that an elastic 
axis exists. For an unswept wing, this 1 b a straight line which undergoes 
no vertical displacement when the wing is subjected to a pure torque par- 
allel to this axis and along which no twisting gradient exists when verti- 
cal loads are applied anywhere along this Line. For an unswept wing of 
conventional construction, this simplification is usually quite accurate. 
For a swept wing an elastic axis may toe defined as a straight line which 
assumes a constant slope over its entire length when a twisting moment is 
applied parallel to this line and which has no twisting gradient when 
vertical loads are applied anywhere along this line. Far aspect ratios 
greater than 5 or 6 and for conventional wing construction, a line can toe 
found on the structure which satisfies this definition reasonably well 
except near the root. It is not uncommon to find an equivalent elastic 
axis at about the 35 or 40 percent chord, a line located aft of the leading 
edge a distance equal to 35 or ^-0 percent of the local chard. 

The assumption of an elastic axis involves the tacit assumption that 
chordwise bending of the wing is negligible. It follows, then, that the 
motion of the wing at any spanwise coordinate can toe described toy two 
coordinates, the vertical displacement of some point on the chord, and the 
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angle of twist of the chord. If wing motion la described In terms of 
vertical motion of the elastic axis and twisting motion about this line, 
then these motions are not coupled through the action of elastic forces 
in the wing except In the root region for a swept wing. 

The root region of a swept wing Is necessarily a relatively compli- 
cated structure. However, for aeroelastic problems an equivalent simple 
structure .can be found which Is completely satisfactory for wings of large 
aspect ratio . This can be demonstrated by the following reasoning. The 
outer sections of a wing exhibit definite beam! Ike properties, but In the 
region of the root considerable warping of the wing surface must take 
place. The aerodynamic forces near the root of the wing are therefore 
not adequately described by strip theory. In addition, the Inertia 
effects of this section are not readily computed. However, the effects 
of the aerodynamic farces on the root section are insignificant for ordi- 
nary flutter computations. This has been demonstrated many times with 
the analog computer by removing the aerodynamic forces on the Inboard 
cell of the finite-difference structure. The inertia forces are also 
insignificant compared with the elastic forces transmitted by the root 
section, and it Is therefore possible to replace this section for purposes 
of analysis by a set of "influence coefficients" relating transmitted 
forces to relative displacement of an outer section of the wing relative 
to the fuselage. It been found that In sane cases these Influence 
coefficients resemble coefficients for a simple beam extending straight 
Into the fuselage and attaching there In sane simple way. The wing 
structural axis then consists of a short section which may be perpendic- 
ular to the fuselage center line and which Is simply attached to a swept- 
back elastic axis which extends to the wing tip. 

Methods for determining the equivalent structure are outside the 
scope of this report. Since this structure varies greatly with the par- 
ticular wing construction vised, it was necessary to choose a simple though 
typical root structure far this study. That chosen is Illustrated In 
figure 3 where the elastic axes are shown by dotted lines. The break In 
the elastic axis Is assumed to be at the edge of the fuselage, and the 
mriH inside the fuselage is assumed to be straight and perpendicular to 
the airp lane center line . The wing Is assumed to be pinned at the side 
of the fuselage. Consequently, all twisting moment la removed at this 
point It Is not necessary to make any assumptions regarding twisting 
rigidity Inside the fuselage. Bending rigidity Inside the fuselage Is, 
however. I m por t ant for symmetric motion. During the past 6 years, exten- 
sive flutter computations have been made with the electric-analog computer 
for commercial and military aircraft as well as for wind-tunnel models 
Including those described In references 3 and 4. In all cases Investi- 
gated, it has been found that relatively large variations in root condi- 
tions have a negligible effect on the flutter characteristics (in the 
sense described In the Introduction). Observed changes In damping were 
usually in the range 0 < | Ag| < 0 . 05 , which has very small effect on 
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flutter speed unless the curve of g against velocity* Is very flat, near 
zero values of g. Needless to say, both symmetric and antisymmetric 
motion of the airplane must be permitted since the flutter characteristics 
for idle two types of motion may be quite different. 

Fuselage stiffness and Inertia properties usually have such values 
that an assumption of a rigid fuselage for symnetrlc motion alters the 
flutter characteristics little. For fighter planes, the error introduced 
Is negligible. For large bombers, the change in flutter speed may be 
appreciable, but It does not alter the trends to be observed upon vari- 
ation of wing properties. It has therefore been assumed In this study 
that the airplane fuselage Is rigid. Tail-surface flexibility does not 
significantly affect wing flutter problems. A rigid tail surface with 
sufficient area to provide satisfactory static stability has therefore 
been assumed. 


Aerodynamic Forces 

For all the flutter computations given In this report, the aero- 
dynamic forces have been simplified by two important assumptions: 

(1) The air flow Is incompressible. 

( 2 ) If the airfoil is divided Into strips perpendicular to the 
elastic axis, then the forces on each strip can be computed as a function 
of the normal component of the alrstream velocity and the motion of that 
strip Independently of the motion of adjacent strips. 

The first assumption Is not required by analog methods in general, 
but Its use greatly increases the rapidity with which data can be obtained. 
Since the purpose of the study Is not to obtain specific accurate flutter 
speeds but to study trends In flutter characteristics, this assumption 
does not seem unreasonable. With regard to the use of strip theory, two 
assumptions are often found In the literature. In using the "alrstream 
method" the wing Is divided Into strips parallel to the alrstream, and 
the forces and moments on each strip are computed as though the wing were 
not swept and the air flow about the section were a two-dimensional incom- 
pressible flow. The aerodynamic coefficients may be taken to be the same 
as those for an unswept wing or may be modified by a factor cos A. In 
applying -the "normal-component method," the wing Is divided Into strips 
perpendicular to the elastic axis. The aerodynamic forces and moments 
are computed as though the effective air velocity were the normal compo- 
nent v cos A, and the forces depend only on the motion of the Individual 
strip and not upon the motion of adjacent strips (except that some small, 
terms may be Included which are proportional to the twisting gradient and 
therefore dependent upon the motion of the nearest strips). A critical 
discussion of the two alternatives is given In reference 5- Ms refer- 
ence recommends use of the normal -component method. 



10 


NACA aw 3780 


Before adopting the second as sumption , an effort was made to find 
some correlation with experimental results. Reference 3 conta ins experi- 
mental flutter speeds for a wind-tunnel model wing with sweepback angle 
equal to about 35°. ThiB angle is sufficient to give an appreciable 
difference in results obtained with the various assumptions mentioned 
above. Ihe section entitled "Finite-Difference Errors" in the present 
report contains the results of computations which show that the normal- 
component method gives results which are as satisfactory as those given 
by any other method used. 

Equations for det ermining the aerodynamic forces by this method are 
given In reference 5- In the equations given there , several terms are 
found whose theoretical Justification is not well established. Obese 
terms (grouped in special brackets on p. 16 of ref. 3) were found to have 
negligible effect on sample flutter computations . It seems reasonable, 
therefore, to omit these terms from computations involved in the present 
trend study. With these emissions and with obvious changes to conform 
to the symbols notation used in the present report, the equations are: 


P * P x + P 2 + P3 

® Mj^ + Mg + 


p ~ - 

P-L =» “2rt(q n )( 2b ) c f^) ^ + (a + 8 tan A) + ^|(a + 0 tan A) 




Pj = -itph 2 ^h - Xqu) 


P^ * -&c(q n )(2b) 
M 1 “ - X 1 P 1 



ci==-iU ^ 0 + 




* 2 t ) + - V) 


tan A 


Mg * — xjP g 
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“5 - + £t+^« -(£ + x^tanA 

The terms are grouped In the order shown for convenience in establishing 
analog circuits. The last tern in is not found in the corresponding 

equation of reference 3* This term is removed (mathematically) by 
insertion of an equal hut opposite term in Pj^ and a similar term in M^. 

It is added to P^ here because the circuits which generate the term 

a + 0 tan A also provide the term j Kg/v^ (& + § tan A), the last part 

of which is not found in reference 5- As is pointed out below, this term 
has a negligible effect so that its inclusion is of no importance, but it 
is indicated in the expression for P-j_ far the sake of conpleteness . 

It should be emphasized that the dynamic pressure q n is based on v n , 
where v n is the velocity component normal to the elastic axis. The 

coordinates a and 0 are both measured in elastic-axis coordinates. 

The symbolism C(bp/v n ) Is used to represent the Theodors en or Wagner 

function. A short discussion of the interpretation of this symbolic 
representation can be found in reference 6. 

All terms found shove can be represented by simple analog circuits 
with the exception of Pj^ and M^. Examination of equations 6 and 7 of 

reference 3 shows that each term in P^ and is similar to (if not 

equal to) a term found in the special brackets. Since the latter terms 
have been omitted, there seems to be no logical reason for retaining 

and Mtj. Inasmuch as their inclusion greatly complicates the analog cir- 
cuits, these terms were also omitted. 

In addition to the finite-difference approximations and those con- 
tained in the assumptions of incompressible flow and strip theory, three 
other aerodynamic approximations should be mentioned. The first of these 
is the failure to modify aerodynamic forces at the wing tip. The delay 
in the growth of lift forces as described by the Wagner or Theodors en 
functions far two-dimensional flow cannot apply near the tip. Indeed, 
both the delay in lift and the magnitude of the lift must go to zero at 
the tip. The extent of the error Introduced depends upon the importance 
of tip forces In flutter computations. Insofar as their location Is 
concerned, these forces are quite Important, but, because of wing taper, 
the magnitude of the total force per unit length diminishes near the tip. 
Since wings of considerable taper are Involved in this Investigation, It 
Is to he expected that the error will he relatively small. The second 
approximation is failure to compute aerodynamic forces properly at the 
root of a swept wing. As stated earlier, the error Introduced by this 
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app r oxima tion Is negligible, since the aerodynamic force for a large 
section of the wing root can he emitted, entirely without an appreciable 
change In flutter speed. The third approximation Is Introduced by the 
necessity of computing the Wagner function (or the The odors en function) 
electrically. This function Is computed using networks shown In refer- 
ence 1 with an error no greater than 2 percent over the frequency range 
or time interval of Interest. 


FIHITE-DIEFHIMCE EffiDBS 


Finite-Difference Structures 

Ho practical methods have been devised for representing general 
continuous structures with continuous electrical systems. The electric- 
analog computer utilizes lumped electrical elements which can. In prin- 
ciple, be used only to construct analogs for lumped mechani cal systems. 
However, as pointed out In references 1 and 2 , it Is possible to repre- 
sent the dynamic characteristics of beamlike structures by a lumped 
structure based upon finite-difference approximations to partial differ- 
ential equations. It Is convenient to call this lumped system a finite- 
difference structure, whether it Is a mechanical model or an electrical 
analog. These references outline the process by which Inertia and stiff- 
ness properties and aerodynamic forces are averaged or replaced by single 
concentrated Inertias, springs, or forces In the finite-difference 
structure. 

It should, be remarked, at this point, that the assumption of a finite- 
difference structure Insures a finite number of flutter roots or exponen- 
tial functions In the transient response, whereas the continuous structure 
has. In principle, an Infinite number. Since the higher frequency roots 
have high damping, it Is only the lower frequency roots that are of 
Interest. Two or three of these may, however, show essentially zero 
damping simultaneously at a given velocity, and it Is sometimes necessary 
to determine the characteristics of several flutter roots. There is 
obviously a lower limit to the number of cells that must be used to obt ain 
satisfactory accuracy, since each cell adds roughly two roots to the 
system. 

There Is little Information In the literature which pertains to the 
accuracy with which such structures represent the continuous system. 
Reference 7 gives data for static-deflection and no rm al-mode character- 
istics of certain finite-difference structures but no Information about 
accuracy of flutter computations. It Is the purpose of this section to 
summarize work at the Analysis Laboratory of the California Institute of 
Technology which was carried out to determine finite-difference errors 
In flutter computations for several specific structures. 
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Using equations for aerodynamic forces based, on two-dimensional 
strip theory and. linear Incompressible fluid flow, several "exact 
solutions" have been obtained far flutter problems. Some of these are 
found In references 8 and 9* These solutions are exact In the sense that 
no further physical or mathematical simplifications are Involved and the 
only errors are Introduced by round-off errors in evaluating transcen- 
dental functions and Infinite series. Solution of these same problems 
by use of finite-difference approximations to partial differential equa- 
tions provides the most practical way of estimating finite-difference 
errors for other configurations for which exact solutions are not obtain- 
able. It is true that, in all cases mentioned above, the airfoil has 
been assumed to have uniform spanwise properties and that in most prac- 
tical cases the airfoil has a significant taper. On the other hand, 
reference 7 contains a study of the finite-difference errors in the 
deflection characteristics and normal-mode properties of both uniform 
and tapered beams. This study showed no unusual differences in these 
properties, and so it is assumed that the results obtained fdr flutter 
of uniform airfoils are typical, of results that would be obtained far 
flutter of tapered airfoils. 

Although much of the work reported in this section was not done in 
the preBent investigation, it is included here since most of it does not 
appear in any readily available publication. 


Uniform Airfoil With Pinned Ends 

A uniform beam with pinned ends will support only sinusoidal modes 
in both bending and torsion. Flutter modes are also of sinusoidal shape 
and it is therefore possible to reduce the flutter problem to an eigen- 
value problem which can be solved with a high degree of numerical accuracy. 
The finite-difference analogs for a plnned-plnned beam likewise will 
support only sinusoidal modes . It is possible therefore to get exact 
solutions for the finite-difference approximations to the continuous 
airfoil. 

The airfoil chosen for this analysis is described in table I. Far 
the continuous wing, the flutter speed and frequency were found to be 
v f = 692 miles per hour and f f = 12. 72 cycles per second, respectively. 

Analysis of the finite-difference structure was carried out using eight-, 
four-, and two-cell divisions between the pinned ends. Results are given 
in table U and figure 4 . For this particular case it is necessary to 
use more than four cells if flutter speed is to be obtained with error 
less than 2 percent. By use of symmetry conditions at the center of the 
beam, it is necessary to use only half this number of cells with an 
electric analog computer. Thus, use of two analog cells gives a theoret- 
ical error of ah cut 2.2 percent, and four analog cells would give an 
error of only 0.6 percent. 
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Uniform Cantilever Wing With Concentrated Maas 

Analytical determination of the flutter speed of a cantilever wing 
is much more difficult than that for a beam with piqned ends. However, 
other investigators have obtained accurate numerical solutions for a few 
configurations. The most Important of these is described in reference 9* 
This case is of Importance for two reasons: It involves several spanwlse 

positions of a large eccentric concentrated mass which has a great effect 
upon the flutter speed; and, for some positions, at least two completely 
different flutter roots can be found. 

Table HI presents the physical characteristics of the airfoil 
analyzed in reference 9* In this reference, the flutter speed and flutter 
frequency were computed for seven mass locations, data for which are 
reproduced in table 17. Since the location of a concentrated mass may be 
Im p o rtant in flutter analysis, and since all points on a finite-difference 
beam are not equally suitable as an attachment point for a concentrated 
mass, it was believed that a comparison of the above data with finite- 
difference solutions was quite important. Unfortunately, similar accurate 
solutions far a finite-difference structure are not readily obtained, so 
it was necessary to use the electric-analog computer to obtain these 
solutions. The resulting comparison therefore contains both finite- 
difference and analog-computer errors. Previous work has Indicated that 
the latter are probably not greater than 1 percent if the Theodors en 
function is represented accurately. 

In this analysis, two slightly different beam analogs were used. 

In both, the heamllke properties were represented by & system of levers 
(transformers), but in one group the lumped forces were applied at the 
junctions of the levers and in the second group the forces were applied 
at the midpoints of the levers. The analog of the Becond group was once 
thought to give a better approximation since it resembles the Bussell beam 
analog discussed in reference 7 . Recent investigation has shown that this 
belief is without foundation, and the Becond analog is now preferred only 
as a matter of convenience for sweptback wings since it provides the wing 
slope directly at the force stations where it is needed far computation 
of aerodynamic forces. In both cases the cantilever condition at the 
root was provided by a half cell at the root, and the forces nearest the 
tip were applied one half cell from the tip. Thus the first group involved 
an integral number of cells, and the second group involved a half Integral 
(integer plus one-half) number of cells. Five cases were investigated; 

2 , 2 ^, 4, 5^, and 6 cells. Since it was shown in the previous section 
that less than 4 cells was of no interest for present purposes, only the 
results of 4, 5jp and 6 cells are presented in this report. 



naca Gar 3780 


35 


la view of the simplicity of the flutter curves shown In reference 9* 
It was expected that data would he taken at only a few spamri.se mass loca- 
tions. However, It was soon found that the flutter characteristics were 
muc h more complicated than anticipated, and data were taken at 24 mass 
locations In the 6-cell case. The flutter characteristics of the wing 
with variable location of the concentrated mass are sketched In fig- 
ure 3(a). As the concentrated mass Is moved outward from the root, the 
flutter speed drops slightly. At a distance about l6 percent of the total 
span from the root a mini mum is reached, and beyond the 25 -percent posi- 
tion the flutter speed rises very rapidly. At the 30-percent position 
the flutter speed for this root has become equal to the flutter speed of 
a completely different root. The flutter speed for this second root drops 
with Increasing spamrise position of the mass making It impossible to 
det ermine with the analog computer the speed for the original root beyond 
the 30-percent position. The flutter speed for the second root reaches a 
minimum with the mass at the 45-percent position, then rises to a very 
high value as the mass Is moved toward the 75 "Percent position. A flutter 
root which Is probably the second Is observed for mass positions near the 
tip, the lowest flutter speed occurring with mass at the tip. It was 
also observed that divergence of the wing occurred whenever the flutter 
speed exceeded about 5*000 Inches per second. Because of divergence. It 
was not possible to measure with accuracy flutter speeds which exceeded 
divergence speed by more than about 50 percent. As a result, flutter 
speeds with mass near the 75-percent span could not be measured. 

The flutter characteristics for the 4-, 5^ ~* and 6-cell structures 

are shown In table V and figure 5(h)- Data for the seven positions ana- 
lyzed In reference 9 are also plotted In the figure. Inspection of these 
curves shows that many more accurate numerical solutions are required, to 
determine the finite-difference errors for all mass positions. In spite 
of the Inadequate numerical data, an attempt was made to draw a smooth 
curve through the known points taken from reference 9. In doing this the 

5^- and 6-cell analog data were used as a guide In determining the shape 

of the curve. This curve, shown In figure 5(a)* has already been dis- 
cussed. It Is realized that a significant error of as much as 2 or 3 per- 
cent may exist In this curve for seme mass positions, but there was no 
other method for obtaining estimated errors for the finite-difference 
structures. With the understanding that the comparison data may be in 
error In some regions, figure 6 was prepared showing the percentage error 
In flutter speed for the various analogs as functions of the mass location. 

For 5^- and 6-cell structures the average errors are about 2 percent. 

It can be readily seen that, although a 4-cell analog gives very satis- 
factory results for the bare wing (mass position 0), it is necessary to 
use more than 4 cells If errors less than 5 percent are required at other 
mass locations . A further discussion of this Investigation will be found 
In reference 10. 
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As a result of this analysis. It vas decided that all flutter com- 
putations made In this trend study would he made using cells to 
represent one-half of the airplane wing. 


Experimental Correlation 

Wind-tunnel tests hare been made of many model structures. It Is 
difficult, however, to find unclassified data In which the structure Is 
completely and accurately described. In the course of this Investigation, 
two cases were found In which a correlation between experimental and 
computed characteristics could be attempted. The first of these 1 b the 
uniform unswept cantilever wing discussed In the preceding section. The 
flutter speed and frequency observed In a wind tunnel are reported In 
reference 9 and a companion report, reference 11. These data are sunma- 
rlzed In table IV, which also contains the computed values of reference 9* 
A better understanding of the correlation Is obtained If the experimental 
data are plotted with the assumed analytic solution. Figure 3(a) shows 
such a comparison. The correlation for thlB case seems unusually good. 


Flutter Speed of a Swept-Wing Model 

Reference 3 gives results of wind-tunnel tests to determine the 
flutter speed of a model wing with sweepback angle equal to 34. 5° • This 
wing had two concentrated masses attached at approximately the 30- and 
80-percent span positions. In an effort to compare the alrstream and 
normal-component aerodynamics for flutter computations, an electrical 
analog was constructed for this wing. For any sweepback angle, it Is to 
be expected that the two methods will give flutter speeds differing by a 

factor of approximately (cos A) 1 / 2 , unless the aerodynamic coefficients 
are modified by the factor cos A In the alrstream method. In which case 
the two methods should give similar results . The principal difficulty 
encountered was determination of the properties of the concentrated masses 
on the wing, since reference 3 does not give complete Information about 
these masses and their geometrical location. The best data that could be 
deduced from, this report are given In table VI. Since the masses are 
alined with the alrstream but are represented In elastic-axis coordinates, 
a product of Inertia between roll and pitch exists. Since no such infor- 
mation was available, the product of inertia was omitted from computations, 
and the rolling Inertia about a chord line was assumed to be one-half as 
large as the pitching Inertia about the elastic axis. It Is believed 
that these approximations and simplifications will affect the results by 
less than 1 percent. 

A comparison of observed and computed characteristics Is given In 
table VH. Qhe first three normal-mode frequencies show satisfactory 
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agreement, with differences of 1, 3/ and 3 percent, respectively. She 
flutter speed computed with either representation of aerodynamic forces 
is lower than the wind-tunnel value. In the case of the airs tr earn method, 
the discrepancy is 19 percent, or, if the aerodynamic coefficients are 
modified, 11 percent. Using the normal-component method, the discrepancy 
is 12 percent. Flutter frequency is in error shout 20 percent in all 
cases. Although the observed differences are relatively large in all 
cases, it is concluded that the normal-component method recommended in 
reference 3 la satisfactory for this model. 


CHARACTEIISTICS OF FOUR REPRESENTATIVE AIRCRAFT 


Plan forms and stiffness and inertia data were chosen after surveying 
the various fighter, bomber, and transport planes developed in recent 
years. Four representative airplanes were chosen, two fighters and two 
large bombers . Smaller attach bombers and transports wqire not included . 
because of lack of time. The airplanes chosen are riot similar in al.l , 
respects to any particular set of four airplanes, but they do have stiff- 
ness and inertia properties which resemble four specific aircraft. Plan 
form, sweepback angle, elastic-axis location, and concentrated-mass loca- 
tions were, however, chosen more arbitrarily so that .this report could 
remain unclassified. Hie four basic plan forms are shown in figure 3* 

The basic fighter A has a bare unswept wing with span of about 500 inches, 
taper ratio of 2.0, and aspect ratio 6. The basic fighter.- Bias a wing 
sweepback angle of 30°, a span of about 400 inches, and a taper ratio 
of 2.0. The two basic wings have the same length measured along the 
elastic axis and the same chords measured perpendicular to the elastic 
axis. 


The basic bomber A has an unswept wing with span of about 
1,700 Inches, taper ratio of 2-5, and aspect ratio 12. It has a concen- 
trated mass representing an engine nacelle at the 0.46-span position 
with center of mass about one-half chord forward of the elastic axis. 

The basic bomber B has a wing sweepback angle of 30°, a span of about 
1,300 Inches, and a taper ratio of 2.4. It also has a concentrated mass 
representing an engine nacelle at the same relative position as for 
bomber A. The two basic wings have the same length measured along the 
elastic axi s and approximately equal chords when measured parallel to 
the air stream. 

Mass per unit length, pitch inertia per unit length, bending rigidity, 
and torsional rigidity were drawn as smooth curves approx! mating the 
characteristics of some typical modern aircraft. As described in refer- 
ence 2, these data must be collected or lumped over distances corresponding 
to the cell leng th of the analog finite-difference structure. The assumed 
curves and the lumped values are shown in figure 7. The lufaped values 
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are also listed In tables VUI to XI, ■which give all pertinent character- 
istics of the basic airplanes. 

Eight Important parameters of the basic airplane wings were varied 
in an effort to find similar features In the flutter characteristics of 
the various wings . The quantities varied and the extent of their vari- 
ation Is sunmarized as follows: 


Quantity varied 


Minimum value Maximum value 


Wing mass density, per unit basic . . 
Wing pitch inertia, per unit basic 
Bending rigidity, per unit basic . . . 
Torsional rigidity, per unit basic . . 
Center-of-mass location, percent chord 
Elastic-axis location, percent chord . 

Chord, per unit basic 

Sweepback angle, deg 


0.5 

2.0 

0.5 

2.0 

0.67 

1.5 

0.67 

1.5 

25 

60 

30 

50 

0.67 

1.5 

0 

1*5 


With the exception of sweepback angle, these quantities were varied 
one at a time from their basic value. However, far an four basic air- 
planes, seme or all of the parameters were varied for two or three values 
of sweepback angle. It is realized that the above variations do not 
constitute a comprehensive survey. However, to a considerable extent the 
changes In flutter speed due to several variations are additive If the 
variations are small and are made simultaneously. Another limitation is 
that the flutter characteristics are affected by the spauwlse variation 
In the first seven quantities listed. The two fighters and two bombers 
constitute four changes in the spauwlse variation of these quantities 
but unfortunately are cases in which four or five of them are varied 
simultaneously. Other quantities which were thought to have second-order 
effects were not considered. Among these are altitude (represented by 
ratio of air density to wing mass), fuselage mass and pitching Inertia, 
and tall configuration. This does not Imply that flutter velocity la 
in dependent of altitude, but with very minor variations the flutter 
velocity varies inversely as the square root of the air density. Sea- 
level air density was used throughout this study. 


It la improbable that b embers of the plhn form and size studied will 
he flown without eng in es on the wing. Consequently, the basic cases of 
interest are those in which a concentrated mass is located there. On the 
other hand, it is of Berne Interest to compare the characteristics of the 
hare wing as well as those of a wing with concentrated mass. Both bomb- 
ers A and B were studied with bare wing as well aa with concentrated 
mass In the basic position on the wing. 


For purpose of reference. It Is necessary to assign a number to 
designate each particular case. The group discussed above comprises 
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173 cases. The assignment of case numbers Is shown in table XU. This 
table shows most rapidly the various cases that were studied. 

Concentrated masses on fighter wings usually consist of fuel tanks, 
bombs, or similar stores. It is Impossible, therefore, to select a 
single value for mass and Inertia which can be regarded as typical. For 
certain positions, many values far mass and inertia were chosen, although 
in most cases the number of values was restricted by the time available 
far computations . For reference purposes, the basic mass far fighter 
planes was arbitrarily chosen to be one-quarter of the mass of the entire 
wing (half of the mass of one side), the pitching radius of gyration was 
set equal to 30 Inches, and the roll radius of gyration was assumed to . 
be 1 5 Inches or less. Specific data for the two fight erB are listed in 
table XHI. 

Concentrated masses for bomber airplanes are usually engine nacelles, 
with a mass which can be predicted within a factor of 2. Nevertheless, 
it is of some Interest to study the effect of various mass values in 
these cases also. Basic mass value far both bombers was assumed to be 
13 pound-seconds squared per inch, which corresponds to a weight of nearly 
6,000 pounds. Pit ching radius of gyration was assumed to be 35 inches. 
Basic mass position was assumed to be at the 0.46-span position and 
60 inches in front of the elastic axis. These data are also tabulated 
in table XIU. 

The concentrated-mass characteristics varied in this study are: 

! l) Mass 

2) Pitching Inertia about carter of mass 

3) Spanwise location 

4) Ghordwlse location 

5) Pitching flexibility 

The assignment of case numbers is more difficult for this phase of 
the study. Although specific spanwise positions were chosen, it was not 
possible to choose chordwlse positions beforehand. The chardwlse posi- 
tions were chosen as the data were obtained. In some cases more than 
20 positions were used for a given spanwise location. Consequently, one 
case number was assigned to all chordwlse variations at a given spanwise 
location. A summary of all variations with the corresponding case num- 
bers is given in table XIV. 

Pit ching flexibility of the concentrated mass was varied in six cases 
involving both bombers. In all cases, the chordwlse location of the cen- 
ter of mass was basic (table XIII ). In three cases the mass was in basic 
spanwise position and in three cases the mass was at the tip. Case numbers 
are given in table XIV. 
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TRENDS IN FIUTTSt CHAEACTHUSHCS 


Reference Quantities and Graphical Presentation 

Results of the study of trends in flutter characteristics -which are 
listed in table XV are given in miles per hour and in per unit values of 
a reference speed. The reference velocity chosen is 

v 0 «= UO^in./sec = 568 mph 

Obviously a flutter speed of 1-5 would not represent a realistic value 
since this would correspond to supersonic speed with a Mach number of 
about 1 . 1 . However, such* a number still has useful significance for 
two reasons: (l) A major purpose of this study is to establish trends 

and to determine what configurations tend to be more or less stable than 
others, and (2) a change in stiffness is equivalent to a change in veloc- 
ity, so that a structure with one-half the stiffness of another, but other- 
wise unchanged, would exhibit a flutter speed l/f2 times as great as 
that of the other, a value equal to 1.06 or 600 per hour in the case 

given above. 

AH geometrical, structural, and inertia quantities are given in 
per unit values. For example, distances are measured in units of the 
airplane semi span and masses, in terms of a basic value. For conversion 
to specific mechanical units, the reference quantities will be found in 
figure 3 j which shows the plan farms, figure 7, which gives inertia per 
unit length and rigidity data for the wings, and tables VTH to XI, which 
list all pertinent characteristics of the four basic airplanes. The 
density of air at sea level was used throughout these computations . The 
value chosen is: 


P » 1. U*6 (lXT 7 ) lb-sec 2 in."^ 

In presenting results graphically, flutter speeds have, in general, 
been reduced to dimensionless values by using as the velocity unit the 
flutter speed of the basic configuration. For example, when plotting 
antisymmetric flutter speed as a function of wing mass density for a 
particular wing such as that of fighter B with A » the flutter 
speeds have been divided by the antisymmetric flutter speed of fighter B, 
A ■= 45 °, with basic wing mass. Bjynraetric and antisymmetric results are 
both presented, rather than choosing the one which gives lowest flutter 
speed. Where such results are presented in the same figure, symmetric 
results are generally indicated by solid lines, and antisymmetric results, 
by dotted lines. Specific numerical values for the flutter speeds and 
flutter frequencies are found in table XV. 
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Mass and Inertia Variations 

In most practical configurations, the normal, mode of vibration with 
lowest frequency is predominantly a bending mode and is usually called 
the first wing bending mode. In the absence of a large concentrated mass 
on the wing, a predominant torsional motion is usually observed in the 
third or fourth mode. Simple flutter can often be predicted with engi- 
neering accuracy using only these two modes as the normal coordinates of 
the structure. When a large concentrated mass is involved, the situation 
is much more complex. Two or more torsion modes as well as two or more 
bending modes become important in flutter computations, and several flut- 
ter roots may be observed which predominantly Involve various ones of 
these modes. Tor eccentric masses it becomes. In fact, impossible to 
speak of bending and torsion modes since many modes will Involve both 
large bending and torsion displacements. 

In those cases in which flutter Involves a bending mode and a higher 
frequency torsion mode, it can be said that a structural change which 
separates the frequencies of these modes ordinarily raises the flutter 
speed, and a change which makes the frequencies more nearly equal lowers 
the flutter speed. It will be observed below that this generalization 
is not always valid. A change in mass density without change in pitching 
inertia has greatest effect on first bending frequency even in cases with 
large sweepback. Consequently, Increase in wing mass density would be 
expected to give an Increase in flutter speed and decrease in mass den- 
sity, a decrease in flutter speed. Changes in pitching inertia would 
normally be expected to have an opposite effect. Such variations were 
made for three fighter configurations, four b are-wing bomber configura- 
tions, and five bomber configurations with concentrated mass. The mass 
density and pitching Inertia were separately changed by factors of 2.0 
and 0.3, making a total of 48 configurations in addition to the 12 basic 
cases. Reference case numbers are given in table XU. 

Tabulation of flutter speed and frequency for each case will be 
found in table XV. The results are also shown In figure 8. As mentioned 
earlier, the flutter speeds have been reduced to dimension less values by 
using as the velocity unit the flutter speed of the basic wing for each 
basic configuration. The trends predicted Shove are found in most cases. 
In the case of b are-wing fighters the effect is very systematic . The 
average of all caseB is given as follows: 

Ifasa-ratae Pitch Inertia mu. Arerage chmge In Tf , 

percent 

0.5 1.0 -12 

2.0 1.0 6 

1.0 .5 8 

1.0 2.0 -9 

The effect, though uniform, is quite small. 
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The results for 'bombers show much less consistency. For cases both 
■with and -without concentrated masses, the effect of wing -mass density 
variation is unpredictable. Nearly half of the cases show trends which 
are opposite to that predicted shove. The addition of a concentrated 
mass at the 0.46-spaa position reversed the trend in several cases. On 
the other hand, change in wing pitching inertia did show a systematic 
trend for all bomber cases . On the average, a change in pitching inertia 
by a factor of two changed the flutter speed about 7 percent. 

The following conclusions can be drawn: 

(1) A change of wing pitching inertia shows a systematic trend 

for all wings, although the effect is small. 

(2) A change of wing mass shows a definite trend far typical 

fighters, although the effect is small. 

(3) Change of wing mass for -typical large bombers with or 

-without concentrated masses shows no systematic trend. 


Stiffness Variations 

It has been pointed out (e.g. , ref. 1, p. 783) that when incompress- 
ible fluid flow is assumed, a change of stiffness is equivalent to a 
change of velocity Insofar as transient response of an airfoil is con- 
cerned. Consequently, it can be said that a uniform increase in stiffness 
will raise the flutter speed by the square root of the factor by which 
stiffness is increased. In most airplanes, it is found that the increase 
in torsional rigidity is primarily responsible far the Increase in flutter 
speed and that, in general, a change in bending rigidity over rather wide 
limits does not change the flutter speed significantly. 

As shown in table XU, 12 configurations were studied to support 
this conclusion. Since both bending rigidity and torsional rigidity were 
separately changed by factors of O.67 and I.50, there are a total of 
I18 case numbers assigned to this group. The results of this study are 
listed in table XV and presented graphically in figure 9* For ease of 
comparison, flutter speeds are converted to dimensionless values, and 
flutter characteristics for changes in bending and torsional rigidity are 
plotted side by side. In general, it was found that change in torsional 
rigidity by a factor of 3/2 or 2/3 increased or decreased the flutter speed 
by 20 percent and that a similar change in bending rigidity had a negli- 
gible effect -upon the flutter speed. Among the 12 configurations studied, 
the following exceptions to this trend were noted: 

(1) Bomber A, A *= 0°: In the antisymmetric case both bending and 

torsional rigidity had roughly equal effects, flutter speed chimgiTig 
±10 percent for the rigidity change given above. 

(2) Bomber B, A = k 5 °i Symmetric case same as case (l) above. 
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(3) Bomber A, concentrated mass at 0.46 span, A » 0°: In the 

antisymmetric case, torsional rigidity* had a 30 percent greater effect 
( +30-percent change In flutter speed) and bending rigidity* had a nega- 
tive effect (710-percent change In flutter speed). 

(4) Bomber A, concentrated mass, A = 30°: In the symmetric case, 

the trend was normal only for Increase in torsional rigidity and decrease 
In bending rigidity*. 

Obese exceptions do not constitute a major deviation, and the trend 
Is considered veil established. 


Local Stiffness Variations 

It Is not to be expected that the same effect will be observed If 
torsional rigidity* Is changed at various stations along the wing. In 
the absence of a concentrated tip mass, any effect on flutter speed must 
vanish for stations near the tip, and presumably the largest effect will 
be observed far stations near the fuselage. Because of the great ease 
with which these data could be obtained, the effect of local stiffness 
variation was obtained far several configurations. 

The analog computer requires lamping or averaging of Inertia and 
stiffness properties. Consequently, it Is possible to determine readily 
only the effect of a stiffness variation which must be assumed to extend 
over the entire length of a cell In the finite-difference structure. The 
basic data consist therefore of step curves. To obtain an approximate 
value for the per unit change In flutter speed per unit change In stiffness 
per unit length at any point along the wing it Is necessary to draw a 
smooth curve passing through this curve such that the areas under the 
two curves are ap p roximately equal. It is believed more suitable to 
present the step curve and let the reader do any smoothing his application 
requires. The configurations studied are listed below: 

(1) Fighter A, A ** 0°, symmetric and antisymmetric, case 10 

(2) Fighter B, A = 0°, symmetric and antisymmetric, case 32 

(3) Bomber A, A * 0°, bare wing, symmetric, case 95 

(4) Bomber A, A « 0°, concentrated mass at 0.46 span, symmetric, 

case 67 

Results of this study are presented In figure 10. Hie abscissa of a 
curve Is the spanwlse station at which the bending or torsional rigidity 
variation Is made. The ordinate is the per unit change in flutter speed 
per unit change in stiffness per unit length along the wing. If, for 
example, the stiffness Is Increased J percent over a distance w along 
a wing of semispan l between the stations d - (w/2) and d + (w/2), 
then the ordinate of the, smoothed curve at the abscissa d/Z when multi- 
plied by j(w/Z) will give the approximate percent change In flutter speed. 
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Data for the two fighters shew great similarity hath for symmetric 
and antisymmetric conditions . The greatest effect Is obtained by changing 
torsional rigidity near the midspan position, slightly outboard for 
fighter A and slightly Inboard for fighter B. She effect of a bending- 
rigidity change was found to be small at all stations. In most cases a 
small negative effect was observed, the flutter speed dropping slightly 
as the bending rigidlly was increased. 

Bomber A without concentrated mass showed a similar trend with the 
following exceptions: ■ 

(1) Maximum Improvement was obtained by changing torsional rigidity 
near the root of the wing ( 0 . 23 -span position). 

(2) Increase In bending rigidity was observed to decrease the flutter 
speed by an amount which was 3 to 10 times greater than that for fighter A. 

Addition of a concentrated mass at sp Bitwise station 0.46 has a great 
effect an this characteristic. The mass chosen Is typical for an air- 
plane engine, and is sufficiently large so that the wing Is to a cert ain 
extent pinned at this point for the particular flutter root Involved. 
Consequently, stiffness changes Inboard of the engine have a negligible 
effect, and changes outboard have an effect very similar to that observed 
for a bare wing of reduced length. 

It should be remarked at this point that the result discussed above 
Is not to be regarded as a trend for all configurations. When the flutter 
Is primarily an outer-wing bending -torsion flutter then this result is to 
be expected. Experience has shown, however, that occasionally an Inner- 
panel torsion mode Is Involved In flutter, and change In torBlonal rigid- 
ity outboard of the nacelle has no significant effect. It Is unfortunate 
that such a configuration was not investigated for this report. 


Center-of-Mass Location 

The location of the wing center of mass has a great effect upon 
flutter speed of an airplane wing. The general trend, Is that flutter 
speed Increases as the center of mass moves forward. It Is not gener- 
ally true that the center of masB is at a constant chord location at all 
spanwlse stations. However, for purposes of studying the trends. It Is 
necessary to assume some basic position for the center of mass. Fast 
experience has shown that a center -of -mass location near the elastic axis 
(usually slightly aft) is both realistic and typical. For this reason 
the basic position of the center of mass was assumed to be the elastic 
axis or 40 percent chord. Variation In center-of-mass location was 
between the 23- and 60 -percent -chord points. Thirteen configurations 
were studied, the various center-of-mass locations comprising 33 cases 
listed In table XU. 
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Hie results are listed In table XV and are shown graphically In 
dimensionless form In figure 11. Pie general trend Is that the flutter 
speed Increases as the center of mass moves forward and decreases as the 
center of mass moves aft, except far center -of -mass locations far behind 
the elastic axis. For positions near the elastic axis, the flutter speed 
changes about 3 percent for a shift In center of mass equal to 1 percent 
of the chord. For the extreme aft positions (60 percent chord) most of 
the curves become quite flat, and In about four cases the flutter speed 
has started to rise slightly as the center of mass Is moved farther aft. 
On the other hand, the curves become very steep for center -of -mass loca- 
tions forward of the elastic axis . In most cases the Increase In flutter 
speed was so great that data could not be obtained for the 23-* and 
32 . 5-pereent-chord locations because the flutter speed greatly exceeded 
the divergence speed. Pie average percentage change In flutter speed 
for a shift In center of mass equal to 1 percent of the chord depends 
upon locatioh of the center of mass as Indicated below: 


Center-of-mass location, percent chord 40 50 60 

Change In flutter speed, percent 3*1 1*7 0.8 


One unusual case was noted. Pie results for fighter B, A *» k-5°, in 
figure ll(a) show an unusual behavior for aft center-of-mass location 
In the antisymmetric case. A study of the frequency of oscillation for 
each position tends to support the conclusion that two different flutter 
roots are Involved. In any case the results are anomalous and could 
bear further Investigation. 

Fighter A, A = 0°, shows another unusual characteristic In the 
antisymmetric case. One flutter root disappears as the center of mass 
Is moved forward of the lf6-percent-chard location. Pile result, shown 
In figure 11(a), Is more easily understood by reference to figure 12 
where the curves of g against v are plotted for this configuration. 

A similar case shown In figure ll(b) has two readily observable flutter 
roots, one with low and the other with high flutter frequency. Data 
for both cases are given In figure ll(b). It Is true that only the one 
with lower flutter speed Is of practical Interest, but for purposes of 
studying trends both are equally Important. A plot of g against v 
for this case Is also shown In figure 12. 


Elastic-Axis Location 

A main component of the aerodynamic pressures on an airfoil Is 
equivalent to a force applied at the quarter chord. Consequently, the 
elastic-axis location relative to the quarter chord determines the 
nature of the coupling between aerodynamic forces and the structure. 

If elastic-axis location alone were changed, both center of pressure and 
center of mass would change with respect to the assumed structural axis. 
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In order to separate the effects due to these two changes, the center of 
mass was moved with the elastic axis In the configurations discussed here. 
Aerodynamic couplin g In which the center of pressure (quarter chord) Is 
forward of the elastic may have a destabilizing Influence while 
center of pressure aft of the elastic axis generally has a stabilizing 
effect. 


Elastic-axis locations between the 30 and 50 percent chord were 
used in the 3 j 6 cases listed In table XU. Results are given in table XV 
and figure 13 - In all cases the expected trend was observed. For an 
elastic axis near the 40 percent chord, the flutter speed changed 
3*2 percent an the average for a shift in elastic axis equal to 1 percent 
of the chord. This effect Is not linear over a wide range, however; the 
flutter speed increases more rapidly as the quarter chord is approached 
and decreases lnore slowly as the elastic axis Is moved aft. For an 
elastic axis at the 30 percent chord, the corresponding change In flutter 
speed was only 1.8 percent. 


Chord Variations 

A change In chard of a wing is usually accompanied by significant 
changes In mass. Inertia, and stiffness as well as changes In other 
characteristics. In an effort to assess the effect of aerodynamic pres- 
sures alone, variations were made in which mass. Inertia, and stiffness 
were held constant while the chord length was changed. Location of the 
elastic axis was maintained at a constant per unit chord station so that 
the distance between quarter chord and elastic axis changed In proportion 
to the change in the chord length. Since the magnitude of the aerodynamic 
force increases with chord length and since the predominantly destabilizing 
lag of the Theodorsen function Increases with chord length. It Is to be 
expected that the flutter speed will decrease as the chord length Is 
Increased. 

Four configurations were studied In which the chord length was 
changed by factors of O.67 and I.30. The eight cases and the configura- 
tions are listed In table XU. Flutter characteristics are given in 
table XV and figure Ik. The results are remarkably uniform. On the 
average, a 7-perceut change in flutter speed results from a 10-percent 
change in chord, smaller chords giving a higher flutter speed. 


Sweepback 

The effect of sweepback upon flutter speed depends upon many factors. 
In conventional wing design, the root structure varies greatly with sweep- 
back angle, and the equivalent elastic axis may show considerable varia- 
tion In position. For wings of large sweepback angle and low aspect ratio. 
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the concept of an elastic axis may- not he useful In describing structural 
properties . From another point of view the problem Is erven more perplexing 
since there Is not general agreement about the nature of the aerodynamic 
forces on a swept wing. In the section of this report entitled "Finite- 
Difference Errors," the results of three methods of computation were com- 
pared with results of wind -tunnel tests of a model wing which was swept- 
back 34. 5°. Two methods were found to give similar results, which were 
significantly better than those of the third. Although the agreement was 
not entirely satisfactory. It was decided to use the aerodynamic forces 
rec ommende d In reference 3* For the present Investigation, the following 
assumptions were therefore made: 

(1) Aerodynamic forces are as discussed In the section entitled 
"Finite-Difference Errors." 

(2) To achieve a sweepback angle, the wing is rotated about a verti- 
cal axis through the Intersection of the unswept elastic axis and the 
side of the fuselage. The tip Is, however, terminated parallel to the 
airstream, so that only the span measured along the elastic axis Is 
unchanged In length. 

(3) Structural properties of the wing are unchanged by sweepback. 

(4) The center of mass of the fuselage is moved aft as the sweepback 
angle Is increased so that It coincides roughly with the center of pres- 
sure of the wing. 

(3) Ho modifications were made for aerodynamic forces at the tip. 

The five basic configurations are shown In table XU, which gives 
reference numbers for the 17 cases. The results are given in table XV 
and figure 13 • Flutter characteristics of the two fighters show a 
reasonable correlation, and. In general, a decrease In flutter speed for 
sweepback angles other than zero. However, the b embers do not show a 
correlation with the fighters or with each other. It Is significant that 
a substantial change In flutter speed with sweepback angle was observed. 

In one case, flutter speed Increased more than 60 percent for a 4-3° sweep- 
back, while other cases showed a 30-percent decrease for sweepback angle 
of about 25°. 

In addition to the cases above. It is possible to crossplot the 
variation of flutter speed with sweepback angle for the following param- 
eter variations of b ember B: Wing mass density, wing pitching Inertia, 

bending rigidity, torsional rigidity, and center-of-mass location. Host 
of these are plotted In figure 1 6. It Is Interesting to note that the 
general trend for bomber B is to a great extent independent of these 
variations . 
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Concentrated-Mass Pitching Flexibility 


The an ginAH on present-day bombers are sometimes mounted In nacelles 
on pylons seme distance below the wing. Because of the inherent flexi- 
bility In such a structure and its fastening to the wing, the dynamic 
characteristics of the ^ gin* are altered. Because of the synmetry of 
the structure. It Is possible to write two sets of equations for the 
nacelle, one Involving pitching, vertical, and fore and aft motion and 
the other involving lateral, rolling, and yawing motion. These sets are 
uncoupled except through interactions with the wing. The characteristics 
represented by the equations involving pitch have a greater effect on 
flutter characteristics, or, stated in another way, the assumption of a 
rigid pylon for lateral motion has not ordinarily been observed to intro- 
duce great differences in flutter characteristics . This assumption 
becomes less valid for wings with large sweepback. On the other hand, 
a significant variation may be observed as the pitching f lexib ilities are 
varied. For pit ching motion it is usually quite accurate to assume an 
effective center of rotation at seme point in front of and below the 
elastic axis. Unless a specific case is to be considered, however, it 
is just as satisfactory to assume this center of rotation at the elastic 
axis, since variation in the location of this point has only a second- 
order effect. Consequently, in this study the center of rotation for 
pitching motion was established at the elastic mein and the pylon was 
assumed rigid far lateral motion. 


Six cases shown in table XU were investigated: 


1 1) Bomber A, 
2) Bomber A, 
3) Bomber B, 
4} Bomber B, 
5) Bomber B, 
6 ) Bomber B, 


A » 0°, mass at 0.46 span, case 176 
A = 0°, mass at tip, case 179 
A «* 0°, mass at 0.46 span, case 177 
A = 0°, mass at tip, case l80 
A = 30°, mass at 0.46 span, case 178 
A » 30°, mass at tip, case l 8 l 


In all of these cases the chordwlse position of the mass was basic, 
60 inches forward of the elastic axis. 


In presenting the results, an e ffort has been made to put the data 
in dimensionless form. Thus, the flutter speed is given as a per unit 
value of the flutter speed with rigid connection. This basic flutter 
speed can be found in table XIV. The flexibility is conveniently meas- 
ured by the normal-mode vibration frequency of the nacelle with the wing 
held rigid in pitch. However, instead of using the value of frequency 
in cycles per second, this frequency is measured in per unit value of 
the flutter frequency with rigid connection. Values for the flutter fre- 
quency with rigid connection can also be found in table XIV. There are 
two frequencies of the nacelle which might be regarded as significant. 

One of these is the cantilever frequency in which the wing is held rigid 
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In "both pitch said plunge. However, for large bombers, the wing has such 
great flexibility in vertical bending that greater significance m i ght be 
attached to the frequency when pitching motion is constrained and vertical 
motion is completely unrestrained. Because of the location chosen for 
the basic mass, the difference in these frequencies is a factor of 2 , the 
frequency with vertical motion unrestrained being higher. For presenta- 
tion of data, this higher value of frequency was chosen, because in those 
cases where a "tuning" effect was observed the maximum effect occurred 
when this frequency was equal to the flutter frequency for the basic 
rigid mass. One exception to this is observed in the discussion below. 

Results are plotted in figure 17. Nine of the twelve cases show a 
predominant decrease in flutter speed as the rigidity is reduced from 
an infini te value. Seven of these cases show a minimum flutter speed 
when the nacelle frequency is nearly equal to the rigid flutter frequency. 
tmh decrease varies between 7 and 37 percent with an average value of 
lB percent. Two cases show a decrease in flutter speed, but no tuning 
effect. The maximum rate of decrease occurs, in fact, when the nacelle 
frequency is far below the rigid flutter frequency. In both cases the 
flutter speed drops to an asymptotic value about six-tenths of the basic 
value. 


Three of the twelve cases show an increase in flutter speed as the 
rigidity is reduced from an infinite value. In two cases increase takes 
place in the region where nacelle frequency is roughly equal to the flut- 
ter frequency. In both cases the flutter speed increases more than 
30 percent. The last anomalous case showB a resonance or tuning effect. 

It is anomalous for two reasons: (l) The flutter speed rises to a sharp 

peak about 10 percent above basic value, and (2) this occurs when nacelle 
frequency is twice as great as the flutter frequency. It should be pointed 
out that, for this rigidity, the flutter frequency is equal to the nacelle 
frequency with wing attachment constrained in bending as well as pitch. 


Effect of a Concentrated Mass 

Many aircraft structures have engines, stores, or external fuel 
tanks attached to the wing in such a way that they act dynamically as 
concentrated masses . It has long been known that the location of such 
a mass has a significant effect on flutter. Unfortunately, other aero- 
dynamic and' seme structural problems do not permi t location of such a 
mass so that maximum flutter speed is obtained. On the other hand, wi thin 
the restrictions Imposed by other considerations, it is often possible 
to Improve flutter characteristics significantly by proper choice of mass 
location. 

This investigation has included a detailed examination of the effect 
of a concentrated mass on the flutter characteristics of several config- 
urations of the four basic airplane wings. Preliminary study of this 
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effect showed such Interesting and unusual effects that the scope of the 
investigation was expanded beyond that originally proposed. The resulting 
data are so voluminous that It Is difficult to present them effectively. 

In particular. It Is Impractical to construct a table which gives all of 
the data obtained, and so graphical presentation Is required. Two methods 
have been adopted In this report. For a given spanwise location of the 
mass, the flutter speed can be plotted as a function of the cbordvise 
location. This has been done for all cases Investigated. Since the con- 
centrated mass Is itHnwi with the alrstream. It Is most convenient In 
cases with sweepback to move the mass parallel to the alrstream rather 
than perpendicular to the elastic axis. 'Where sufficient data are avail- 
able, these curves can be summarized In a single diagram In which line s 
of constant flutter speed are shown on a drawing of the wing plan farm. 

Far the concentrated mass located anywhere on such a contour line, the 
flutter speed will be the same. Die result Is essentially a topographic 
map of the fliitter -speed surface, Tiber e each point on the plan form 
represents a possible location for the concentrated mass. 

Several difficulties arise with both methods of presentation. The 
main source of difficulty -lies In the fact that several important flutter 
roots exist far a wing with concentrated masB . For certain locations of 
the mass, one root will show lowest flutter speed, while for other loca- 
tions another root will have the lowest flutter speed. Since the analog 
computer Is essentially an electrical model. It Is usually impossible to 
find one of these flutter speeds if another root has a flutter speed far 
below the first. It is possible, therefore, to find with certainty only 
those portions of a given flutter-root surface which lie beneath all 
other flutter-root surfaces. Far one configuration studied, four such 
distinct surfaces were positively Identified and It was not possible to 
establish that surfaces appearing at widely separated regions were or 
were not related. In most cases the roots were differentiated by obtaining 
essentially marginal stability far two distinct roots along the line 
where the two surfaces Intersect. It can be readily appreciated that many 
points are required to establish the flutter-speed contours, particularly 
where several intersecting surfaces are Involved. It was. In fact, impos- 
sible in the time available to obtain sufficient data to establish all 
Interesting features about these contours. However, It Is believed that 
all important features are shown correctly In the figures presented here. 

Hie curves which show flutter speed as a function of chordwlse 
position at a fixed span frequently show Intersections between different 
flutter roots. In identifying these roots, it is useful to know the 
flutter frequency associated with each root. Hie simplest way to present 
these data Is to show the value of frequency at a few selected points 
along each curve. Where roots Intersect and both frequencies were meas- 
ured, both values are shown. In some cases where actual frequencies ware 
not measured, low, medium, or high frequency are shown. 
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For convenience, the concentrated mass was always placed at the 
center of a finite-difference cell excerpt In three cases where additional 
Information vas obtained by placing It halfway between cells. It is 
convenient to Identify these locations by the cell number as has been 
done in table XIV, which assigns a case number to each configuration. 

It must be remembered, however, that the cell divisions are slightly dif- 
ferent for bomber and fighter airplanes, and therefore the spanwlse sta- 
tion for a given cell number will be different. The location of these 
stations In terms of unit span Is given In table XIH. In the figures, 
the spanwlse position Is correctly given as a fraction or per unit value 
of the wing semispan. 

The size of the concentrated mass and its pitching and rolling 
Inertia also affect the flutter speed. Since past experience has shown 
that rolling Inertia has a small effect, a few cases were chosen for 
further Investigation of the magnitude of this effect. For a concentrated 
mass located In the wing it Is reasonable to assume a radius of gyration 
which Is a small fraction of the average half chord. For a mass suspended 
below the wing, it Is unlikely that the distance will exceed half of the 
average half chord. Two- values for radius of gyration were chosen, equal 
to 0.1 and 0.5 times the average chord for the fighter planen. In all 
cases considered there was no significant difference In flutter character- 
isticfl when the rolling inertia was varied frcm zero to the maximum value. 
The variation was. In fact, so insignificant that none of the data is 
presented In this report. In what follows it may be Assumed that the 
rolling Inertia of the concentrated mass has any value between the above 
limits. Since the mass of the concentrated mass and its pitching Inertia 
have a greater effect. It la necessary to consider variations In these 
quantities In several typical cases. The basic values for mass and 
pitching Inertia (|or radius of gyration k) have been discussed in the 
section entitled ''Characteristics of Four Representative Aircraft" and 
are given In table XTU. The variations of these values are summarized 
In table XIV. 

The most logical way to give the results Is to present first the 
flutter characteristics for the basic mass an each particular airplane. 
Five airplanes were chosen: 

! l) Fighter A, A ■» 0° 

2) Fighter A, A - 45° 

3) Bomber A, A = 0° 

4) Bomber B, A ** 0° 

5) Bomber B, A = 30° 

Since It was difficult to choose a typical radius of gyration for a mass 
on a fighter plane, three values were used. These values are 6, 15, 
and 30 Inches, as shown In table XIV. Figures 18(a) and 18(b) show the 
effect of chordwlse location of the mass at five spanwlse positions far 
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fighter A vith A * 0 °. Results for all three radii of gyration In pitch 
are plotted on the same sheet using different symbols for each value. 
Circles are used for the smallest value, k = 6 Inches} triangles are 
used for k * 30 Inches; a solid line with no symbols Is used far the 
Intermediate value. Abscissas for all curves are chordwlse distance 
from the elastic axis measured as per unit value of the wing semispan. 
Similar data for fighter A, A = 4 - 5 °, are shown In flgureB 18 (c) and 
18(d); six spanwlse stations were used In this case. One surprising 
feature can he noted In all of these figures: The characteristics are 

relatively Independent of the pitching inertia, even though this Inertia 
Is varied from a very large value (k = 30 Inches) to nearly zero 
(k = 6 Inches). This does not mean that at any particular point the 
flutter speeds are identical, hut the overall shapes of the curves show 
remarkable similarity. 

Although these figures (figs. 18(a) to l8(k) ) give a good picture 
of the flutter characteristics. It is easier to Interpret the results If 
all data are combined to construct flutter contours as discussed earlier. 
Such contours for the minimum and mm-timim values of k are shown In 
figures 19 (a) to 19(h). These figures support the following conclusions: 

(1) A chordwlse position aft of the elastic axis Is almost always 
undesirable. 

( 2 ) The 30 - to 30 -percent-span position and the tip location are 
generally undesirable. 

( 3 ) A position forward of the elastic axis and near the 70- to 
80-percent-apan position will. In general, greatly Increase the flutter 
speed. 

Since time did not permit a complete study of the characteristics 
far fighter B, data were obtained only for the cases of 0° and 4-3° sweep- 
hack with mass at the tip. Comparison of the results shown In fig- 
ure 18(e) with the corresponding data far fighter A in figures 18 (a), 
18(h), 18(c), and 18(d) shows that for this location there does not seem 
to he any significant difference In characteristics. Whether It is safe 
to extrapolate this result to other mass locations cannot he said at 
this time. 

Since the study of fighter A shewed that the pitch radius of gyration 
had a small, effect, and since the pitch radius of gyration of a bomber 
engine Is relatively well defined, it was decided to use only one value 
In the study of bomber airplanes. However, the practice of using one 
engine or two engines on a single pylon, as well as the different sizes 
of engines, gives a possible variation In mass which might well exceed 
a factor of two. All bomber data were therefore obtained with both basic 
mass and half basic mass. For ease of comparison of the two sets of 
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data, tjaey are plotted side by side In the figures. Data for bomber A , 

A ■= 0°, are contained In figures l 8 (f ) and 18(g). Results for b ember B, 

A 3 o° are given In figures 18(h) and l 8 (i), and the case of bomber B, 

A =■ 30°, Is summarized In figures l 8 (j) and l£(k). Again It Is possible 
to simplify Interpretation of these figures by combining the results 
Into flutter-speed contours. However, It can be seen that data for basic 
mass and half basic mass are very similar, and so such contours have been 
prepared only for the cases with basic mass. The flutter-speed contours 
are shown In figures 19(i) to 19(n). A study of these figures showB some 
deviations from the results for fighter airplanes. The following con- 
clusions can be drown: 

(1) A position aft of the elastic axis Is almost always undesirable. 

( 2 ) The tip region Is generally undesirable as a location for the 

mass. 


(3) With few exceptions, any position forward of the elastic axis 
and between the root and the 90 ^percent-span position will give flutter 
speed equal to or greater than the bare -wing flutter speed. 

(4) However, there are. In most cases, no practical locations which 
give any significant Improvement In flutter characteristics . Two cases 
will he noted In which the speed might he Increased 40 percent. The 
others are restricted to a 10- or 20-percent Improvement. 

Since fighter planes showed remarkable variation In flutter charac- 
teristics with mass position, It was believed necessary to examine the 
effect of changes In the size (mass) of the concentrated mass. This was 
first studied at two spanwlse positions far fighter A, A *= 0°. The 
positions are the tip and station 3 (O .79 span). Flutter characteristics 
as functions of chordwlse position were measured for several values of 
mass. The three values of radius of gyration given In table XIII were 
used far all mass values, except that for very small values of mass only 
the 6 -Inch value was used. However, because of the similarity of results, 
data are presented only for the maximum value (k = 30 inches) and minimum 
value (k =* 6 Inches). Case numbers are listed in table XIV. 

The results for tip location shown In figure 20 (a) show a very 
interesting, progression In characteristics as the mass Is reduced to 
zero. Most striking Is the fact that no significant change takes place 
when the mass is varied frcm twice basic value to half basic value. Sven 
with mass reduced to 8 percent of basic value, the three flutter roots 
far antisymmetric motion and the two flutter roots for symmetric motion 
can still he Identified, though their characteristics are by this time 
somewhat altered. Similar data far the mass at 0-79 span are presented 
In figure 20 (b). Figure 20(a) concluded Illustrates graphically the 
danger in extrapolating results. For a mass at the 0.10 chordwlse 
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position and symmetric motion, a mass equal to 0.008 basic mass increases 
flutter speed to I.07. Doubling the mass Increases it to l.l6. Again 
doubling it vill increase the speed to same unknown value greater than 
1.40. However, if the mass is again increased by a factor of about 2.8, 
the flutter speed (of another root) vill have dropped to 1.00 again. 

Flutter characteristics were also measured for basic and half basic 
mass at the wing tip with sweepback angle of 45° • This case was chosen 
because the antisymmetric characteristic for a 30-inch radius of gyration 
(fig. 18(d)) showed a very unusual characteristic. Figure 2l(a) shows 
that reduction, of the mass by a factor of two eliminates the anomalous 
behavior, but in all other respects gives results which are essentially 
the same as the basic mass. Figure 2l(b) gives similar data for the 
mass at the O.79 span position with A - 45° • Again the results for 
basic mass and half basic mass are not significantly different. It is 
perhaps unwise to attempt any general statement, but there is every indi- 
cation that the essential features shown in the flutter-speed contours 
of figure 19 would not be profoundly altered if either mass or pitching 
inertia were increased or decreased by a factor as large as 2.0. 


Accuracy and Probable Errors 

A brief statement about the expected accuracy of these results has 
been purposely deferred until the end of this report. It is believed 
that this discussion win be more meaningful after the reader has observed 
the nature and scope of the data obtained. The analog computer 1 b not 
composed of perfect electrical elements. For example, the inductors used 
in this study have loss characteristics corresponding to a damping factor 
of about g * 0.01. Transformers also have significant losses. The 
electrical analog of the airplanes studied in this report had an electri- 
cal damping corresponding to a structural damping between g = 0.02 and 
g - 0.03. Ibis is not greatly different from the damping to be found in 
conventional aircraft construction, so no corrections were made for this 
internal damping. 

No general statement can be made about the effect of random computer 
errors. Some give rise primarily to an error in the damping factor of 
the roots, in which the (g - v) curve is shifted vertically. Other errors 
give rise basically to an error in velocity, in which the curve is prima- 
rily shifted horizontally. Since the slope of the curve of g against v 
is by no means constant, it becomes impossible to give a specific figure 
for accuracy of flutter speed. In unusual cases, where roots are of the 
type shown in figure 12, flutter may be predicted when in fact it will 
not occur far this root at all. Oils distinction, which must be made 
mathematically, is of no importance in practical cases. An airplane which 
shows a flutter damping of g =* 0.01 and is therefore theoretically stable 
is not to be regarded as any more satisfactory or useful than one which 
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shows a damping g = -0.01 and would therefore theoretically fly apart. 
Manufacturing tolerances and the safety factors required In aircraft 
will not permit use of an aircraft uni ess It Is moderately stable for 
a significant variation In all structural parameters. It is believed 
that, exclusive of errors Introduced by the finite-difference structure 
and approximations In the aerodynamic theory, the results obtained In 
this study have a probable error In damping factor of about g ** +0.02 
or a probable e rr o r In flutter speed of about 2 percent, whichever is 
applicable In the light of the above discussion. However, trends obtained 
by variation of structural parameters are considerably more accurate than 
this would Imply, since any error would persist with roughly the same 
value In all cases Involving such parameter variations. 

For cases 1 to l8l, it Is possible to construct curves of g 
against v although they have not been made a part of this report. 

These curves serve to show the steepness with which the roots pass 
through flutter and give sane Indication of the accuracy of the flutter 
speed. This situation does not exist for cases l8l to 289. For 
concentrated-mass variations, computations were carried out In such a 
way that only flutter speeds and frequency were Obtained. Consequently, 
it Is not possible to delineate areas which are "safe (l) 11 from the standpoint 
of flutter. It Is known, for example, that where a long pendant lobe Is 
observed, as for case 235 tn figure lB(d), the system Is barely unstable 
everywhere within this lobe. There ere similar regions elsewhere, for 
example case 212, which Is antisymmetric (fig. 20(b)), where the system 
is barely unstable within an elliptical-shaped boundary. Similarly, 
case 2 lA 1b barely stable In this region and yet no flutter root Is even 
shown since the system, does not become actually unstable at any point. 

These remarks are not made to show the flutter curves to be value- 
less, but to caution the reader against making Inferences not contained 
In the report and not legitimately supported by the data presented here. 


SUMMARY OF RESUIZDS 


Certain trends In flutter characteristics for typical modern air- 
craft seem to he Indicated by this study of the Incompressible flutter 
characteristics of aircraft wings. In some cases a few deviations are 
found. It- Is probable that If more extreme aircraft designs were con- 
sidered, even more would be observed. Nevertheless, this summary may 
serve as a useful guide. 

(l) In the following table are listed the average changes In flutter 
speed for a 1-percent change In each parameter, the change being made 
from the basic value except for center -of -^nas s and elastic-axis location. 

In which cases several locations are assumed. 
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Change in parameter 

Average change in 
flutter speed, percent 

Increase wing mass 1 percent 

0 

Increase wing pitching Inertia 1 percent 

-.1 

Increase bending rigidity 1 percent 

0 

Increase torsional rigidity 1 percent 

•5 

Increase wing chord 1 percent 

-.7 

Center of mass forward 1 percent of chord from - 


40-percent location 

3.0 

50 -percent location 

1-7 

60 -percent location 

.8 

Elastic axis forward 1 percent of chord from - 


ltO-percent location 

3.2 

30 -percent location 

1.8 


( 2 ) Local iz ed change In torsional rigidity is most effective in 
changing flutter speed of a hare-wing airplane if the change is made 
between midspan and root. Tor bombers with large concentrated mass on 
the wing, torsional rigidity either inboard or outboard of the mass will 
usually govern flutter speed, depending on the type of flutter existing. 
The effect of sweepback was not observed to have a systematic effect. 
Pitching flexibility of the concentrated-mass support has a definite 
influence an flutter speed. In many cases a tuning effect was observed, 
with a ID- to 40-percent decrease in flutter speed. This effect was not 
always observed; in some cases, the flutter speed was significantly 
Increased. 

( 3 ) Perhaps the most Interesting results will be found in the effects 
of a concentrated-mass location. Tor a wide range of mass and radius - 

of -gyration values the results were very systematic. Tor fighter-type 
planes it was found that: 

(a) Aft cbordwlse positions are usually undesirable. 

(b) The 30- to 50-percent-span and tip locations are generally 

undesirable. 

(c) A forward location near the 70 - to 80-percent-span position 

will, in general, greatly Increase flutter speed. 

Tor bomber-type planes theBe results are somewhat modified: 

(a) Aft cbordwlse positions are usually undesirable. 

(b) The tip location is generally undesirable. 
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(c) With few exceptions, any location forward of the elastic axis 
and between root and 0.90-span position Is satisfactory, although flutter 
speed Is rarely greater than bare-wlng flutter speed by any significant 
amount. 

There Is perhaps no need to remark that these results can be altered 
by Introduction of a flexibility In the concentrated-mass support. 


California Institute of Technology, 
Pasadena, Calif., June 6, 1955* 
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TABLE I.- PHYSICAL CHARACTERISTICS OF A PIHBED-PIHBED AIRFOIL 


Length, In 288 

Half chord, h. In. 48 

Hass per unit length, m, Ib-sec 2 /sq. In 0.01035 

Inertia per unit length, I, Ib-sec 2 5*178 

Bending rigidity, EE, Ib-in. 2 1.412 (lo9) 

Torsional, rigidity, GJ, lb-in. 2 6.87(10®) 

Xp In 9*8 

x^. In -9.6 

AJLr density, p, Ih-sec^ln.^ 0.0845(l0“^) 


TABLE H.- COMPUTED FEUTTER SPEED AND FREQUENCE QF FINLEE-DIFFEREECE 

PUSHED-PIH11ED AIRFOIL 


Dumber of cells 

E9 

A*f/v» 

f f , 

cps 

A- * 

00 

692 

0 

12.7 

0 

8 

688 

.006 

12.6 

.008 

4 

677 

.022 

12.3 

*032 

2 

644 

.069 

11.0 

.134 
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SABLE IH.- PHX5ICAL CHARACTERISTICS OF UKEFOEM CANTUSVHl WING 


Half chord, b, In. . 

Span, Z, in. 

Mass per unit length, m, li-8ee 2 /sq. In. . . . 
Pitching Inertia per unit length, lb-sec 2 . . . 

Flexural rigidity, SI, (lb)(s^ In.) 

Torsional rigidity, GJ, (lb)(sq. In.) 

Elastic-axis position, Xq, In 

Center of mass position, X5, In 

Mass of concentrated mass, lb-sec 2 /ln. . . 

Pitch Inertia of concentrated mass, lb-sec 2 /in. 
Center of mass of concentrated mass, In. 

Air density, p, lb-sec 2 /in.^ ......... 


4 

48 


1.877 (lO- 4 ) 
8.oo(io J 0 
i.4ot(io5) 
0.692(105) 
. . -O .504 


O.156 


8.23(lO" 5 ) 
. . 0.1636 
. . -3-274 

1.155(10"^) 
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TABLE 17.- THEORETICAL AND EXPERIMENTAL FLUTTER CHARACTERISTICS OF 
CANTILEVER WING WITH CONCENTRATED MASS 
[Data taken from references 9 and 11 1 


Mass 

location, 
per unit 
span 

Calculated 

Experimental 

Vf, 

mph 

per unit 
(a) 

f £> 

cps 

mph 

v f , 

per unit 
(a) 

ff, 

cps 

0 

227 

0.400 

25-27 

228 

0.401 

22.1 

.167 

— 



221 

.388 

19.1 

.229 

226 

•397 

19-23 

221 

.388 

17.4 

.292 




233 

.410 

16.3 

.333 

WMM 



256 

.451 

15.5 

•354 

277 

.488 

28.04 

260 

.458 

* 16 . 3 - 26.8 

.625 

c 359 

.631 

30.68 

... 


(d) 

•938 

273 

.481 

25.67 

261 

.459 

(d) 

•959 

251 

.442 

24.87 

251 

.442 

21.8 

•979 

--- 



231 

.407 

21.6 

1.000 

205 

.360 

23.60 

218 

.384 

21.4 


a 1.0 per unit velocily Is ID** - In. /sec or 568 mph. 

*Tbis experimental record seems to show nearly simultaneous diver- 
gence and flutter at two frequencies. 


Calculated divergence speed Is about 279 mph. However, 
speed can still be calculated mathematically. 

^Divergence was observed experimental ly. 


a flutter 
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TABIE VI.- PHYSICAL CHARACTERISTICS OF SWEPTBACK WINS WITH 

CONCENTRATED MASS 

(Data taken from reference 3* Naas of wing la for portion outboard of 
root restraint. More detailed information will be found in refer- 
ence 3- Data for concentrated masses are not given explicitly in 
reference 3 and must be regarded as only approximate f] 


Wing characteristics: 

a Span, in 48.3 

^Root half chord, bp, in 3*2 

^Tip half chord, b^, in 2.36 

Wing total mass, %, Ib-sec 2 /in 0.00784- 

Tunnel fluid density, p, l±>-sec 2 /in.^ 3-4o(lO"7) 

Sweepback angle. A, deg 34.3 

Concentrated-mass characteristics: Inboard Outboard 

Mass, me, lb-sec 2 /in 0.00806 0.00452 

Fitch inertia about elastic axis, li-sec 2 /in. . . 0.0712 0.0192 

Per unit spanwlse position (from root) 0.30 O.78 

Center of mass position, (x^ c , 121 "1.74 O.50 


fa) Measured along elastic axis. 

(b) Measured perpendicular to elastic axis. 


3jype of frequency 

Hormal mode frequencies, 

ops j art - 

Mode 1 

Mode 2 

Mode 5 


Experimental model frequency 

6.97 

30.9 

37-9 

Measured analog frequency 

6.91 

52.6 

39.1 


Flutter characteristics 

nype of result 

y £> 

mpb 

£sr*hr*> 

*£> 

cps 

Wind-tunnel results, Vf e 

193 

0 

20.1 

Analog results, noarmal-canponant method 

170 

-.12 

24.2 

Analog results, alrstream method 8. 

157 

-19 

24.0 

Analog results, alrstream method* 3 

■ 173 

-.11 

24.0 


a Lift coefficient - 2*. 

^Lift coefficient ■ 2* cob A. 
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TABIE VIH.- CHARACTER1BEECS OF BASIC FlCEffiCER A 
(a) Physical Characteristics 


Sweepback angle. A, deg 0 

Semispan of wing, a l. In 238 

Cell size for finite-difference structure. Ay, in. 34 

Boot chord, 13 bp. In 106 

Up chord,' 0 b x . In 33 

Taper ratio 2.00 

Aspect ratio 6 

Wing elastic axis, percent chord 40 

Wing center of mass, percent chord 40 

Total wing mass external of fuselage, m^, U>-sec 2 /ln 10 . 7 

Fuselage mass, m^, Ib-sec 2 /in 21 

Fuselage radius of gyration, pitch, 0 In 100 

Fuselage radius of gyration, roll, in 23 

Fuselage center of mass aft of elastic axis. In. 0 

Tail center of pressure aft of elastic axis. In 230 

Tall area, sq. In 3,000 

Air density, p, li-sec 2 /ln.^ 1. 146 (lO"?) 


(b) Inertia and stiffness values lumped far finite-difference structure 


Station 

n 

2 

3 

n 

5 

6 

Per unit span, 8, 

a 

Cl 

• 

0 

0.357 

0.300 

0.643 

0.786 

* 

• 

0 

Half chard, 13 b. In. . . . 

47.3 

43.5 

39.8 

36.0 

32.2 

28.4 

Lumped mass, 2, 

lb-sec 2 /in 

1-73 

1.34 

0.98 

0.66 

0.40 

0.23 

lumped pitch Inertia . . . 

695 

568 

405 

253 

•126 

53 


d W 1D J (dy/BI) .... 

19.8 

61.3 

1 12 

205 

415 

m 


d lX) 1 0 y*( dy/GJ) .... 

56.3 

145 

219 

365 

667 


850 


a Measured along elastic axis. 

^Measured perpendicular to elastic axis. 

°About elastic axis. 

^Stiffness values are lunped between mass stations. 
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TABLE IX.- CHARACTERISTICS OF BASIC FESTER B 
(a) Physical characteristics 


Sweepback angle. A, deg 30 

Semispan of ving, a 1, In 238 

Cell size for finite -difference structure, £y, in 34- 

Boot chord, 13 bp. In. 106 

Tip chord, 11 bg. In 33 

Taper ratio 2.00 

Wing elastic axis, percent Chord . 4-0 

Wing center of mass, percent chord 4-0 

Total vlng mass external of fuselage, lb-sec^/in 14.06 

Fuselage mass, mf, Ib-sec 2 /in. ......... 21 


Fuselage radius of gyration, pitch, 0 In. . . 

Fuselage radius of gyration, roll. In 

Fuselage center of mass aft of elastic axis, in. 
Tall center of pressure aft of elastic axis. In. 
Tall area, sq. In 

Air density, p, lb-aec^in.^ 


100 


25 

0 

230 


. . . 3,000 
1.146 GxrO 


(b) Inertia and stiffness values lumped far finite-difference structure 


Station 

fl 

2 

3 

a 

5 

6 

Per unit span, 8, 

0.214- 

0.357 

0.500 

0.643 

0.786 

0.928 

Half chord, 13 b, in. . . . 

lumped mass, m. 

47*3 

43.5 

39.8 

36.0 

32;2 

28.4 

Xb-sec2/in 


1.43 

1.32 

1.15 

0.94 

0.66 

Lumped pitch Inertia . „ . 

j 

577 

503 

411 

310 

204 


d 10*-° J (dy/KL) 

33.0 

100 

156 

221 

321 

447 

— 

d 10l0 J (dy/GJ) 

71.4, 

193 

283 

'551 

400 

436 

224 


s Measured along elastic axis. 

1) Measured perpendicular to elastic axis. 

°About elastic axis. 

^Stiffness values are lumped between mass stations. 
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ZABLE X.- CHARACTERISTICS OF BASIC BOMBER A 
(a) Physical characteristics 


Sweepback angle. A, deg 0 

Semispan of wing, 8, l , In 845 

Cell size for finite-difference structure, Zy, In. 130 

Root chord,* 5 by, In 200 

Tip chord,* 5 b^, In. 80 

Taper ratio 2.5O 

Aspect ratio 12 

Wing elastic axis, percent chord 40 

Wing center of mass, percent chord 40 

Total wing mass external of fuselage, j H-sec 2 /in 39*7 

Fuselage mass, m^, Ib-sec 2 /in. 120 

Fuselage radius of gyration - pitch, 0 In 240 

Fuselage radius of gyration - roll. In 30 

Fuselage center of mass aft of elastic axis. In. 0 

Tall center of pressure aft of elastic axis. In. 700 

Tall area, sq In 20,000 

Air density, p, Ib-sec 2 /ln.^ 1.1)46 (lO“T) 


(b) Inertia and stiffness values lumped for finite-difference structure 


Station 

n 

2 

3 

D 

5 

6 

Per unit span, 8, 

0.154 

0.308 

o.46i 

0 

• 

ON 

0.769 

0 

VD 

Half chord,* 5 b, in. . . . 

90.8 

8I.5 

72.3 

63.1 

53-9 

44.6 

lumped mass, m, 

U>-sec 2 /ln 

6.28 

5.46 

4.16 

2.47 

0.98 

0.50 

Lumped pitch Inertia . . . 

13,300 

10,600 

7,150 

3,640 

1,170 

360 


*10 10 J (dy/El) .... 

3.3 

8.4 

14.3 

30.5 

101 


— 

d lX)10 J (dy/GJ) . . . 


23.1 

30.9 

47-5 

143 

1 

618 


a Measured along elastic axis. 

^Measured perpendicular to elastic axis. 
c About elastic axis. 

Stiffness values are lumped between mass stations. 
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TABLE XI.- CHARACTERISTICS OF BASIC BOMBER B 
(a) Physical characteristics 


Sweepback angle. A, deg 30 

Semispan of tring, a l, in. 845 

Cell size for finite-difference structure, Zy, in. 130 

Root chord, 13 hr, in 170 

Tip chord,*’ b x , in. 70 

Taper ratio 2.4-3 

Wing elastic axis, percent chord 40 

Wing center of mass, percent chord 40 

Total iring mass external of fuselage, j lb-sec 2 /in 48.7 

Fuselage mass, m^ , Ih-sec^/ln. 120 

Fuselage radius of gyration, pitch, 0 in 240 

Fuselage radius of gyration, roll, in. 50 

Fuselage center of mass aft of elastic axis, in. 150 

Tail center of pressure aft of elastic axis, in. 700 

Tall area, sq in 20,000 

Air densl-hy, p, Xb-sec 2 /in-^ 1.146 (iO"T) 


(b) Inertia and stiffness values lumped for finite difference structure 


Station 

a 

2 

3 

a 

5 

6 

Per unit span, 8, 

0.154 

0.308 

0.461 

0.615 

0.769 

0.923 

Half chord, 13 b, in. . . . 

77.3 

69.6 

61.9 

54.2 

46.5 

38.8 

Lumped mass, m, 

lb-sec 2 /in 

6.28 

5-23 

4.26 

3.35 

2.74 

2.47 

Lumped pitch Inertia . . . 

3,940 

3,080 

2,200 

1,560 

910 

390 


a io 10 J~(ay/EI) 


18.8 

31.0 

54.9 

91.6 

141.0 

4 ioio J (ay/oj) 

E 

15.5 

22.6 

42.2 

77.8 

110.0 


a Measured along elastic axis. 

^’Measured perpendicular to elastic axis. 

°About elastic axis. 

^Stiffness values are lumped between mass stations. 
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TABLE Xm.- C0BCEKIBA3ED-MASS CHARACTERISTICS AND IDCATEDNS 

(a) Characteristics 



Fighter 

Bomber 

Basic mass, Xb-sec 2 /in 

2.68 

15.0 

Fitch radius of gyration, 8. In. . . . 

50 

55 


15 

17.5 


6 

7 

Roll radius of gyration. In. . . . 

15 

0 


Basic spanwlse position,* 1 In. . . . 


590 

Basic chordwise position, 0 In. . . . 


60 


a About elastic axis; alrstream coordinates. 

^Outboard from center line measured along elastic axis. 
c Forward of elastic axis, parallel to alrstream. 


(h) Location 


Airplane 

Locations of 

concentrated mass, per unit 

span, a 

at station - 

K9 

2 

5 

D 

Q| 

5 

5.5 

6 

Tip 

Fighter 

0.2114- 

0.358 

0.900 

0.643 


0.786 


0.929 

1.00 



Bomber 

.154 

.308 

.461 

VO 

• 

0.692 

.769 

0.846 

.925 

B9 


a Distances are measured In per unit span along elastic axis. 
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XABX2 XV.- TEtJTTEa SPEED AHD JHK4JJHCT 




Etjmetrlc 

Antjayetrlc 

Case 

Airplane 

Tf» 


ff. 


▼f. 

n 



nph 

per unit 

epa 

■ph 

per unit 

El 1 

1 

Fighter A 

989 

1 - 7 1 * 

lb. 2 

1,091 

1-92 

8.7 

2 


1,060 

1.90 

10.1* 

1,102 

1 . 9 b 

6.13 

5 


892 

1.37 

16.8 

1,068 

1.88 

01.3 

4 


835 

1 A 7 

12. b 

1,0U 

1.78 

8.1 

5 


1,102 

1.9b 

13-0 

l,lb2 

2.01 

8.7 

6 


972 

1-71 

lb. 8 

1,083 

1.91 

9-1 

7 


1,000 

1.76 

13-7 

1,109 

1-97 

8.3 

8 


1,222 

2.15 

16.8 

1,369 

2-bi 

10.2 

9 


795 

1.1*0 

12.0 


886 

1.56 

7 -b 

ll 


813 

1 . 1*3 

lb.o 

.98 9 

1.7b 

8.6 



710 

1-25 

13.9 


636 

1.12 

16.6 

12 



1892 


8-3 

13 


619 

1.09 

12.9 

< 

I 5 b 5 

[778 

.96 

1-37 

16.0 

7-9 

Ob 


773 

1.36 

lb.i 


838 

1.31 

8.2 

W 

lfi 


.793 

1.V 

2.28 

12.0 


790 

1-39 

7.1 


938 

V 4 I 

lb.b 


733 

1.29 

17.2 

17 

l.Gj 

IB 


9^9 

I.67 

15.0 


R 

1.32 

17.7 

19 


1,028 

1.81 

05-7 


1.33 

17.3 

20 


1A15 

2 . 1*9 

l *.5 

1,267 

2.25 

16.1 

21 


8b7 

1 - 1*9 

lb. 9 


78b 

1.38 

17.0 

22 


790 

1-39 

lb. 6 


795 

l.bo 

03-9 

23 

Fighter B 

l,0rft 

1.89 

7 -b 

1,017 

0-79 

6.9 

* 


\,006 

1-77 

5-3 

i, 

$ 

1.06 

b.8 

25 


>p 

1.?* 

10.8 


1.65 

9-2 

26 


909 

1 - 71 * 

6.6 


926 

1.63 

6.1 

27 


1,031 


8.0 

1,07b 

1.89 

7 -b 

28 


1,023 

8.2 

1,017 

1.79 

7.2 

29 


1,01*0 

iL.es 

6.8 

1,017 

0-79 

6.8 

30 


1,273 

2.2b 

8.b 

i,2bb 

2.19 

8.3 

31 


835 

1 . 1*7 

6.7 


830 

l.b6 

5-9 

33 


881 

1-35 

8.0 


898 

1.58 

6.5 

3 b 


767 

1-35 

8.1 


82b 

l-b 3 

6.1 

35 


655 

1-15 

7-3 


TOO 

1.25 


36 


l,Shb 

2.19 

7.8 

1,227 

2.16 

7-3 

51 


852 

1.50 

7 -b 


813 

l.b3 

6.3 

38 


781 * 

1.38 

7-1 


727 

1.28 

3.3 

39 


1.338 

2.39 

8.5 

1,392 

2 .b 3 

8.3 

bo 


9*3 

1.66 

8.1 


R 

1.32 

11.2 

tt. 


932 

1.6b 

9-0 


o-bi 

oa.b 

1|2 


983 

0-73 

9-2 

1,023 

1.80 

11.3 

b 3 


1,031 

1.83 

6-7 

1,031 

1.99 

9.2 

bb 


®J7 

l.b 9 

11.8 


8b 7 

l.b 9 

03.9 

b 3 


861 * 

1-52 

8.0 

1,080 

1.90 

1D.0 

b6 

*7 


■a 

1.89 

0..73 

9.8 

10.2 

1,059 

9 ® 

2.0b 

1-73 

12.6 

12.9 

b8 


960 

1.69 

8.7 

1,07b 

1.69 

11.2 

b 9 


1,176 

2.07 

11.7 

1,303 

2.31 

13.8 

50 


801 

1.1*1 

8.3 


801 

l.bl 

10.6 

31 


2,171 

1,188 

3.82 

11-5 

2,037 

3.62 

13.8 

5 ® 


2.09 

8.8 

l.lb2 

2.01 

12.0 

53 


898 

a? 

9-3 

1,017 

1.79 

11.2 

5 **- 


835 

8.9 

1,307 

2.30 

7-9 

53 


756 

1.33 

8.b 


989 

l-T 1 * 

o.b 

36 


1 , 31*7 

2.37 

10.2 

1,230 

2.20 

02. b 

37 


781 * 

1.38 

8.0 


9 b 9 | 

I.67 

n.b 
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9KBB XT.- HUBS BRED AID IHEglHTT - CanfelznaA 


MB 


tymatrlo 

Antliji—iitrlo 


Alxplana 

r t > 

T f» 


T f» 

T f' 




aph 

par unit 

atm 

nph 

par unit 

asm 

58 

Bariber A 

£ 

3-25 

9-* 

869 

1-33 

8.8 

59 


1.10 

7-5 

6*8 

1.1* 

7-0 

& 


9*5 

1.66 

12.2 

& 

1.71 

lfi.8 

61 


602 

1.06 

8.2 

1.20 

7-3 

62 


£ 

1.58 

10.2 

1,0*0 

1.® 

9.9 



1.21 

9-7 

761 

1.3* 

8.6 

A 


B 

1.29 

1.58 

iii 

1,000 

1,227 

I.76 

2.16 

12.6 

66 


565 

•99 

7.9 

619 

1.09 

7-° 

68 


625 

1.10 

6.8 

676 

3-19 

8.0 

69 


585 

1.05 

8* 

J»3 

.99 

7-* 

70 


5*0 

•» 

7-8 

435 

.80 

6.7 

71 


55* 

•9* 

0.6 

J699 
. 1835 

-95 

7-7 

72 


1,085 

1.91 

17.9 

1.23 

1**7 

B 



1,195 

2.10 

17.7 

f 886 

{l,011 

1.36 

3-78 

10.6 

1B.9 



•78* 

1.58 

8.6 

655 

1.15 

7-5 



1.165 

2.05 

20.9 

,905 

1-39 

20.1 



966 

I.70 

l*-5 

J597 

1TP5 

at 

8.2 

15-* 



1,017 
1 ,03* 

1-79 

21.8 

778 

1.37 

9.9 



1.82 

17-* 

7** 

1.31 

9-7 



866 

3-56 

19-5 

Eg 

1-23 

9.0 



1,085 

1.91. 

25-7 

1-53 

U.0 



8*7 

1 -*9 

1*.2 

f « 

1.07 

7-9 

02 


750 

1.52 

18.5 

]e*2 

I676 

1.13 

3-19 

h 

85 


580 

1.02 

17.5 

1363 

\608 

•99 

1.07 

Vo 

A 


818 

1.** 

9.* 

886 

1.36 

9.9 

85 


86* 

3U52 

7.8 

767 

1.53 

8.7 

86 

a 


£ 

1.5* 

l-l? 

UL* 

8.0 

9® 

86* 

1.73 

1.52 

9-6 

7-1 

88 


909 

1.60 
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Figure 1.- Transient response corresponding to roots of characteristic 

equation. 



Figure 2.- Variation of g with velocity for typical flutter roots. 
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Spanwise position of mass 


Figure 


(a) Theoretical and wind-tunnel flutter characteristics. 

5*- Flutter characteristics of uniform cantilever wing with con- 
centrated mass. 
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(t>) Flutter characteristics of finite-difference analog. 
Figure 5.- Concluded. 
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Figure 6 .- Flutter-speed errors of finite-difference analog of uniform 
cantilever wing with, concentrated mass. 






Per unit span 


(b} Basic fighter B. 
Figure 7 *- Continued, 
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Per unit span 

(d) Basic bomber B. 
Figure 7*- Concluded. 












Bariber A. 

8.- Continued 





Wing mass density 


Wing pitching Inertia 


(c) Bather B. 
Figure 8.- Continued. 
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(c) Bomber Bj 
Tigure 9*- 


Symmetrlc 


Antisymmetric 



wing. 



Per unit torsional rigidity 


(&) B ember 
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Figure 11.- Flutter speed as function of center-of-mass position In per- 
cent chord. 
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Symmetric 

Antisymmetric 




(c) Bomber B. 
Figure 11.- Concluded, 
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a C of M 60% Case 13 


(a) Fighter Aj A = 0°; antisymmetric 


g o 



o C of M 40% Case 72 
O C of M 50% Case 82 
a C of M 60% Case 83 


(h) Bomber A j A = 50° j concentrated mass on wing; antisymmetric. 

Figure 12.- Flat of dam p in g factor g against velocity for seme un u sua l 

flutter roots. 
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Percent chord 


Figure 13.- Flutter speed as a function of elastic-axis position in per 

cent chord. 















Figure lk.~ 


-Symmetric 
- Antisymmetric 



speed os function of chord length. 
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Symmetric 

Antisymmetric 



Figure 15- - Flutter speed as function of erweepback. angle. 
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Figure 16.- Effect of sweepback on ■various modifications of ’b omb er B. 



(a) Mass at 0.^6 span. 

Figure 17 .- Flutter speed, as a function of c cncerrtrated-maa s pitching 

frequency. 
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Mass at tij>. 
17. - Concluded 
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(a) Fighter A; A = 0°; symmetric. 

Figure 18.- Flutter characteristics against chordwlse position of 

concentrated mass. 
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-2 -.1 0 .1 .2 .3 .4 

Chordwise position 


(a) Concluded. 
Figure 18.- Continued. 
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(h) Fighter A; A = 0°; atrtisynmetrlc. 


Figure 18.- Continued. 
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(to) Concluded. 
Figure 18.- Continued. 
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• Figure, 18.- Continued. 
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(d) Continued. 
Figure 18 .- Continued. 
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(g) Bomber A} A = 0°; antisyimietric. HP, high frequency. 
Figure 18.- Continued. 








naca. rar 3780 


99 



(g) Concluded. 
Figure 18.- Continued. 
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(h) Bomber B; A = 0 °; symmetric. 
Figure IB.- Continued. 
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(i) Bomber B; A = 0°; airtisynmetric. 
Figure 18.- Continued. 
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-.1 0 .! .2 -J 0 .1 .2 

Chordwise positron 


(l) Concluded. 
Figure 18 .- Continued. 
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(j) Concluded. 
Figure 18 .- Continued.. 
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(k) Bontoer B; A = 30°; antisyranetric . IF, low frequency. 
Figure 16 .- Continued. 
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(k) Concluded. 


Figure 18.- Concluded. 
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Forward 


(a) Fighter A; A = 0°; symmetric; radius of gyration, 6 Inches; 
1.0v.p = 989 miles per hour. 


mmm 


Forward 




(h) Fighter A; A = 0°; symmetric; radius of gyration, 30 Inches; 

1 . 0 v f = 989 miles per hour. 

Figure 19 .- Contours of constant flutter speed with concentrated mass 









(&) Fighter A; A = 0 °; antisymmetric; radius of gyration, 30 Inc h es; 

l.Ovj * 1*091 miles per hour. 



Figure 19.- Continued. 
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(e) F ight er A; A = 45 °; synmetric; radius of gyration, 6 inches; 
l.OVjf = 1,028 miles per hour. 


Figure 19.- Continued. 
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(a) 'Fighter A; A = i +-5°; antisymmetric; radius of gyration, 6 inches; 

1.0v f = 881 miles per hour. 

Figure 19.- Continued. 
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(n) Bomber B; anti symmetri c } 1.0v f = 1,244 miles per hour. 


Figure 19 .- Concluded. 
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(fc) Concluded. 

‘ Figure 20 .- Continued. 
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Chord wise position 


Chordwlse position 


Figure 21. 


(a) Basic mass and 0.5 Basic mass at tip. 

Comparison of characteristics for basic mass and. 0.5 basic 
mass. Fighter A, A = 45°. 
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